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STUDY  OF  VACUUM  ULTRAVIOLET  GENERATION  BY 
COHERENT  RESONANT  FREQUENCY  MIXING  IN  METAL  VAPORS 

I •  GENERAL  OBJECTIVES 

The  general  long  term  goal  of  this  project  is  to  demonstrate 
an  efficient  method  of  frequency  up-conversion,  using 
simultaneously  "resonant"  and  "coherent"  enhancement.  Conditions 
of  multiphoton  resonances  are  sought  to  enhance  the  nonlinear 
suscept ibi 1 ity  responsible  for  harmonic  generation,  coherent 
propagation  effects  are  exploited  to  minimize  depletion  of  the 
fundamental  by  the  (multiphoton)  resonant  absorption  and  by 
photoionization,  and  to  maximize  the  distance  over  which  phase 
matching  is  maintained.  We  have  demonstrated  experimentally  the 
theoretical  maximum  conversion  efficiency  of  0.5%,  for  the  full 
beam  cross  section,  in  lithium  vapor  (section  IV).  A  theoretical 
study  has  been  made  in  preparation  for  experiments  of  5th 
harmonic  generation  in  condition  of  four  phonon  resonant  coherent 
interaction  in  mercury  vapor  (section  V). 

This  project  leads  also  to  a  analytic  technique  to  measure 
nonlinear  susceptibilities.  The  measurements  of  third  harmonic 
versus  phase  difference  between  a  pair  of  exciting  pulses 
presented  in  section  IV  enabled  us  to  measure  for  the  first  time 
the  phase  angle  of  the  nonlinear  susceptibility. 

Harmonic  generation  is  only  one  of  many  applications  that 
can  be  enhanced  by  multiphoton  resonant  coherent  interactions,  it 
is  shown  in  section  VI  that  a  technique  of  multiple  pulse 
sequences  can  be  used  to  achieve  spectral  selectivity  with 
ultrashort  pulses,  even  in  broadband  molecular  systems. 


II.  BASIC  APPROACH 


With  the  availability  of  shorter  and  shorter  pulses,  hence 
higher  peak  intensities,  one  might  expect  simultaneous 
improvements  in  nonlinear  optical  techniques  such  as  sum  and 
difference  frequency  generation.  However,  as  the  pulse  become 
shorter,  their  increased  bandwidth  will  interact  with  a  larger 
number  of  levels  of  the  atom,  leading  to  very  complex  propagation 
characteristics,  and  increasing  losses  of  energy  through 
multiphoton  absorption.  As  the  pulses  become  shorter  than  the 
phase  relaxation  times  of  the  media  involved,  there  is  a 
possibility  to  reverse  the  energy  transfers  between  matter  and 
radiation,  and  return  to  the  fields  the  energy  onc6  lost  by 
absorption.  The  purpose  of  this  study  is  to  find  a  way  to  take 
advantage  of  the  reversibility  of  light  matter  interaction  with 
ultrashort  coherent  interaction,  to  enhance  the  efficiency  of  sum 
frequency  generation.  The  general  method  sought  is  to  "prepare" 
the  medium  with  a  first  pulse,  and  use  a  second  pulse  to  convert 
back  to  the  radiation  field,  at  the  desired  wavelength,  all  the 
energy  lost  to  the  medium  by  the  first  pulse.  In  this  particular 
project,  the  nonlinear  medium  is  excited  by  a  sequence  of  two 
pu 1 ses  . 


One  may  already  wonder  at  this  point  whether  it  makes  sense 
to  divide  the  available  power  into  small  fragments.  Is  the  total 
excitation  not  just  the  sum  of  that  by  each  of  the  successive 
pulses?  Mot  if  the  pulses  are  shorter  than  the  phase  memory  of 
the  matter  -  collision  time  or  internal  phase  relaxation  time  -. 
The  relative  phase  of  the  pulses  determine  which  way  the  energy 
is  to  flow  -  that  is  whether  there  will  be  absorption  or 
stimulated  emission. 

There  is  a  quantum  jump  in  complexity  in  going  from  one  to 
two  pulse  excitation  -  because  one  is  adding  a  dimen  ion  to  the 
problem:  the  phase.  The  experimental  requirements  are  rather 
stringent:  the  source  has  to  be  able  to  generate  high  power 
tunable  pulses  with  a  uniform  wavefront  -  an  essential  condition 
in  order  to  be  able  to  define  the  phase  -  .  One  needs  to  be  able 
to  adjust  the  relative  phase  and  delay  between  the  two  pulses  of 
the  sequence.  The  phase  or  frequency  modulation  of  the  pulses 
should  be  small,  or  it  would  not  be  possible  to  define  an  average 
phase  for  each  pulse.  Finally,  resonant  coherent  interaction  is 
strongly  dependent  upon  the  pulse  shape  and  frequency,  which  are 
not  particularly  easy  to  determine  with  picosecond  pulses. 
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III.  OVERVIEW  OF  THE  SYSTEMS  UNDER  INVESTIGATION 

The  relevant  energy  levels  of  the  two  atomic  system  that 
have  been  the  object  of  this  investigation  are  sketched  in  Fig.  1 
below. 

Li  Hg 


3u> 


Figure  lb 

Relevant  energy  levels 
of  mercury. 

The  first  and  simplest  resonant  enhancement  to  be  considered 
for  harmonic  generation  is  the  two-photon  resonance.  As  an 
example,  we  have  chosen  to  conduct  experiments  on  the  2s  -  4s 
transition  in  lithium  vapor  (figure  la),  for  which  all  the 
relevant  parameters  have  been  calculated  ■'■.  This  atomic  system  is 
(two-photon)  resonant  for  radiation  at  571  nm.  The  proper  partial 
vapor  of  magnesium  will  phase  match  the  third  harmonic  and 
fundamental.  The  energy  density  required  to  significantly 
populate  the  4s  state  is  of  the  order  of  50  mJ/cm  .  Some  of  the 
experimental  difficulties  to  overcome  are: 

a)  Generation  of  pulses  of  6  psec  duration,  1  mJ  energy, 
with  a  flat  wavefront,  a  Gaussian  beam  profile,  and  no 
temporal  phase  modulation. 

a)  Splitting  of  these  pulses  into  a  sequence  of  two 
identical  pulses  of  well  defined  relative  delay  and 
phase . 

e)  Monitoring  the  phase  difference  between  the  two  pulses  of 
the  sequence. 

b)  Preventing  temporal  and  spatial  phase  distortion  in  the 
amplifier  and  in  the  heat  pipe. 


Figure  la 

Relevant  energy  levels 
of  lithium. 


c)  Production  of  an  homogeneous  mixture  of  Li  and  Mg  near 
1,000  °C. 

d)  Finding  a  "transient  phase  matching  condition"  that  is 
entirely  different  from  the  phase  matching  condition  at 
moderate  intensities  (or  in  the  absence  of  coherent 
interactions ) . 

A  two-photon  resonance  in  a  four  photon  process  -  as  the  one 
involved  in  third  harmonic  generation  -  is  a  very  particular  case 
of  resonant  interaction.  Indeed,  as  shown  in  reference  ■*■,  the 
pulse  energy  is  the  main  parameter  determining  the  relative 
population  transfers.  For  pulses  shorter  than  the  phase 
relaxation  (elastic  collision  time),  the  time  is  only  a  scaling 
parameter  (for  instance,  shorter  pulses  will  produce  the  same 
amount  of  third  harmonic  than  longer  pulses  of  the  same  energy). 
This  is  non  longer  the  case  when  higher  order  resonances  are 
involved . 

The  next  system  under  investigation  is  mercury  vapor  excited 
in  conditions  of  four  photon  resonance  (figure  lb).  There  are  two 
resonances  that  are  accessible  to  amplified  picosecond  dye 
lasers:  the  6s  -  6d  transition  (four  photon  resonance  for 
radiation  at  561  nm),  and  the  6s  -  7s  transition  (four  photon 
resonance  for  radiation  at  625  nm).  The  case  of  a  four-photon 
resonance  is  much  more  complex  than  that  of  a  two- photon 
resonance,  because  Stark  shifts  of  various  order  become  important 
at  the  intensities  required  for  effective  four-photon  excitation. 
In  addition,  the  pulse  duration  becomes  an  important  parameter, 
since,  at  constant  pulse  energy,  the  population  transfers  will  be 
largest  for  the  shortest  pulses. 

In  the  particular  case  of  the  6s  -  6d  transition  of  mercury, 
the  6p  level  is  close  to  3/4  the  energy  of  the  6d  level, 
resulting  in  a  resonant  enhancement  of  the  Stark  shift  of  the  6s 
-  6d  transition.  This  Stark  shift  is  so  important  that,  for 
femtosecond  pulses  of  5  J/cm^  tuned  to  1/4  of  the  6s  -  7s 
transition  frequency,  the  6s  -  6d  transition  becomes  also  near 
resonant  1 


IV.  COHERENT  INTERACTIONS  IN  LITHIUM  VAPOR 
IV. 1.  The  Laser  Oscillator-Amplifier  System 


We  have  constructed  a  laser  osci  llator -a mpl  1  f  i e r  system  * 
capable  of  producing  single  pulses  of  1  raj  energy,  at  a 
repetition  rate  of  10  Hz,  tunable  around  570  nm.  The  amplifier  is 
sketched  in  Fig.  2  below.  Design  details  can  be  found  in  Appendix 
A  (ref.  2).  The  first  two  stages  are  pumped  t ransversally ,  and 
the  last  one  longitudinally.  Coupling  between  the  first  stages  is 
minimized  by  synchronizing  the  first  stage  at  the  leading  edge, 
while  the  second  stage  is  synchronized  at  the  trailing  edge  of 
the  pump  pulse.  The  amplified  spontaneous  emission  (ASE)  after 
the  last  stage  is  approximately  1%  of  the  pulse  energy,  without 
the  use  of  any  other  isolator  than  a  direct  vision  prism  between 
the  two  first  stages.  Several  dyes  and  solvents  have  been 
investigated,  which  have  a  maximum  gain  and  a  minimum  absorption 
at  570  nm.  The  purpose  of  the  small  absorption  for  the  unpumped 
dye  is  to  have  the  amplifier  transparent  to  the  unamplified  beam. 


3-STAGE  AMPLIFIER 


(POL  =i) 


CELL  1  AND  2  AT  BREWSTER  ANGLE 


2  Sketch  of  the  Laser  Amplifier 
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an  essential  condition  to  perform  the  critical  alignments  of  the 
interferometers  used  for  the  analysis  of  time  and  space 
coherence.  The  combination  of  Hexaf luoroisopropanol  (HFI)  with 
Rhodamine  6G  shifts  sufficiently  the  absorption  and  fluorescence 
spectra  of  Rhodamine  6G  to  meet  the  conditions  of  no  absorption 
and  maximum  fluorescence  at  570  nm.  In  view  of  the  high  costs  of 
HFI  mixtures  of  water  and  HFI  were  investigated,  which  were 
reported  to  combine  the  spectral  shift  of  HFI  and  the  thermal 
conductivity  properties  of  water  3.  Unfortunately,  the  desired 
spectral  shift  were  not  found  in  the  water  solutions,  which 
separated  into  two  unmiscible  components  after  a  few  days.  For 
operation  at  570  nm,  the  latest  choice  of  dye  solutions  is: 


1st  stage: 
>r.d 

3rd 


2nd  stage: 


!tage: 


Beam  profile 


Rhodamine  6G  .18  mM/1  in  Ethanol. 

Rhodamine  6G  .14  mM/1  in  Ethanol. 

Rhodamine  6G  in  Water  and  Ammonix  0.01  mM/1. 

Careful  position  of  the  pump  beams  relative  to 
the  amplified  beam  in  the  last  two  stages 
provides  complete  control  over  the  spatial 
intensity  distribution,  from  "square"  (fig. 
3a)  to  Gaussian  (fig.  3b). 


SPATIAL  PROFILE  OF  THE  AMPLIFIES  BEAM 
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Figure  3a 


Beam  Profile 


Figure  3b 


seccr:o  ii/irkonic  signal 


>atial  Coherence  of  the  amplified  pulse  is  demonstrated  by  achieving 

uniform  extinction  of  the  beam  in  a  Mack  Zehnder 
interferometer . 

mporal  Coherence  The  measurement  of  temporal  coherence  is  crucial  to 

the  experiments  of  coherent  propagation,  and  has  been 
used  to  test  the  laser  sources,  the  scanning  techni¬ 
que,  and  the  data  acquisition  system.  Measurements  as 
the  one  shown  in  Fig.  4  indicate  that  the  amplifier 
does  not  introduce  any  phase  modulation.  The  ampli¬ 
fied  pulse  can  best  be  represented  by  a  Gaussian 
temporal  profile  with  a  phase  modulation  <^(t)  propor¬ 
tional  to  the  pulse  intensity  ( <p  =  1.43  E^).  This  is 

» 

the  first  time  that  the  technique  of  interferometric 
autocorrelations  5  has  been  applied  to  the  chirp 
determination  of  synchronously  mode-locked  lasers. 
The  small  phase  modulation  that  we  measure  corres¬ 
ponds  to  an  increase  of  bandwidth  duration  product  of 
only  10  %,  hence  is  too  small  to  be  determined  by 
conventional  spectral  measurements. 


Figure  4 

Interferometric  autocorrelation  of  the  6.2  psec  pulses 
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IV. 2.  The  Lithium-Magnesium  Heat  pipe 

Three  successive  prototypes  had  to  be  constructed  to  evolve 
from  the  basic  heat  pipe  design  of  Vidal  and  cooper  ^  to  an 
instrument  which  confines  a  homogeneous,  uniform  and  transparent 
mixture  of  lithium  and  magnesium  vapors,  and  that  can  be  operated 
continuously  for  several  days  at  1000  °C.  The  rather  slow  data 
acquisition  process  (limited  essentially  by  the  slow  repetition 
rate  of  the  laser  osc i 1  la t or -a mp 1 i f i e r  chain)  imposes  severe 
constraint  not  only  on  the  source  stability,  but  also  on  all 
other  components  of  the  experiment.  Obviously  the  heat'  pipe  has 
to  be  perfectly  sealed,  but  in  addition  the  outgassing  of  the 
walls  has  to  be  continuously  and  selectively  (i.e.  without 
pumping  the  He  buffer  gas)  eliminated  during  the  experiment.  The 
selective  pumping  was  achieved  with  cryogenical ly  cooled 
molecular  sieves.  Another  problem  is  a  condensation  of  solid 
lithium  at  the  edges  of  the  hot  zone.  The  resulting  accumulation 
over  several  hours  of  an  annular  ridge  of  solid  has  been 
sufficient  to  obturate  the  optical  path,  in  the  first  two 
prototypes.  The  problem  of  accumulation  of  lithium  near  the  ends 
is  alleviated  with  the  larger  final  version  of  the  heat  pipe,  but 
it  still  limits  its  continuous  operation  to  roughly  one  week. 

Phase  matched  third  harmonic  generation  measurements  have 
been  made  in  homogeneous  mixtures  of  lithium  and  magnesium.  The 
relatively  low  vapor  pressure  of  solid  magnesium  causes  it  to 
escape  (by  sublimation)  from  the  central  region  of  the  heat  pipe 
during  the  heating  period.  Once  the  lithium  has  melted,  the 
magnesium  is  trapped  with  the  lithium  vapor  in  the  central 
region.  Attempts  to  introduce  the  magnesium  in  the  hot  heat  pipe 
(after  the  lithium  has  melted)  have  been  abandoned  as  too 
hazardous . 

A  fourth  heat  pipe  which  consists  in  an  upper  and  lower  plan 
circular  electrode  with  a  central  mesh,  separated  by  a  pyrex 
cylinder,  has  been  constructed.  This  heat  pipe  is  intended  for 
measurement  of  photoionization. 
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The  Data  Acquisition  System 


The  experimental  set  up  is  sketched  in  figure  5.  The  beam 
from  the  oscillator-amplifier  is  sent  through  a  "pulse  sequencer, 
which  splits  each  pulse  into  two  pulses  of  continuously 
increasing  relative  phase  and  delay.  The  two  pulse  sequence 
energy  is  monitored  by  D  3  before  being  sent  through  the  heat 
pipe.  The  transmitted  signal(s)  is  (are)  monitored  by  D4.  All 
detectors  at  the  fundamental  laser  wavelength  ( ,  D3,  and  d4  are 
combination  of  frequency  doubling  KDP  crystals,  filters,  and  UV 
detectors.  In  addition  to  discriminating  the  amplified  pulses 
from  the  nonamplified  pulse  train  (which  has  a  higher  average 
power),  chis  second  harmonic  detection  scheme  provides  more 
accurate  transmission  measurements. 


A  Mack  Zehnder  interferometer7  is  used  to  create  the 
sequence  of  two  pulses  with  continuously  increasing  phase  and 
delay  (as  needed  for  the  experiments  of  coherent  propagation,  as 
well  as  for  the  analysis  of  pulse  coherence).  The  (collimated) 
beam  produced  by  the  amplifier  is  split  into  two  parts  by  an 
adjustable  beam  splitter  (based  on  frustrated  total  internal 
reflection).  One  beam  is  given  a  fixed  delay  by  a  combination  of 
fixed  prisms  optically  contacted  together,  while  the  other  beam 
is  sent  through  an  adjustable  optical  delay  line,  before  being 
recombined  with  the  first  one  with  another  adjustable  beam 
splitter  . 

4 

The  phase  and  delay  are  really  one  single  parameter.  As  the 
delay  is  scanned  continuously,  the  phase  increments  by  2n  for 
every  delay  increment  of  a  wavelength  of  .57  micron.  With  a  pulse 
repetition  rate  of  10  pps,  the  translation  speed  of  the  delay  arm 
of  the  interferometer  should  be  less  than  .1  microns/s  in  order 
to  have  at  least  20  data  points  per  period!  Such  a  low  speed  is 
incompatible  with  a  smooth  continuous  motion.  Because  of  solid 
friction  in  the  drive,  the  delay  stops  moving,  then  jumps  a 
wavelength  or  more,  making  it  impossible  to  maintain  and 
determine  the  delay  and  phase  with  interferometric  accuracy.  The 
solution  is  to  maintain  the  translation  speed  high  enough  for 
continuous  motion,  to  monitor  continuously  the  phase,  and  to 
store  the  data  points  of  successive  periods  (for  example  100 
periods)  into  an  "equivalent"  period.  The  exact  value  of  the 
phase  is  determined  by  monitoring  the  transmission  of  a  He-Ne 
laser  through  the  same  delay  line.  This  monitoring  is  essential, 
because  even  the  best  synchronous  motor  drive  for  the  delay  line 
does  not  provide  a  continuous  speed  for  the  translation  stage, 
because  of  the  errors  introduced  by  reduction  gears.  The  exact 
optical  delay  X  in  the  pulse  sequencer  is  determined  (with  an 
accuracy  better  than  0.03  micron)  from  the  (cos  )_i  of  the  He-Ne 
laser  intensity  recorded  by  the  detector  D2* 

The  amount  of  averaging  that  can  be  done  is  directly  related 
to  the  duration  and  coherence  of  the  pulses.  The  criteria  are 
those  of  the  slowly  varying  approximation:  the  averaging  can  only 
be  performed  over  the  number  of  cycles  over  which  the  pulse 
amplitude  and  phase  can  be  considered  to  be  constant.  With  our 
6.2  psec  pulses,  the  maximum  number  of  periods  over  which  the 
data  could  be  averaged  is  at  most  N  =  300.  As  the  delay  is  being 
continuously  scanned,  the  data  are  accumulated  for  a  total  delay 


increment  of  Xmax  =  i;X  (where  X  is  the  wavelength  of  the  pulses), 
and  sorted  according  to  their  relative  phase: 

=  2  7T  (  X  -  nX)  /x 

The  accuracy  of  the  averaging  process  is  directly  related  to  the 

accuracy  in  the  determination  of  the  wavelength.  If,  for 

instance,  we  choose  to  divide  the  2rr  phase  interval  in  2  5 

segments,  and  average  the  data  over  N  =100  periods,  the  maximum 

allowable  error  in  the  wavelength  determination  is  X/2500.  Fig.  6 

below  shows  an  example  of  such  an  averaging  over  100  periods,  of 

the  second  harmonic  of  the  signal  created  by  the  pulse  sequencer, 
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as  recorded  by  D3,  normalized  to  the  second  harmonic  of  the  pulse 
from  the  laser  source,  as  recorded  by  .  In  that  particular 
recording,  the  data  points  were  sorted  in  25  phase  sectors. 


Figure  6  Recording  of  signal  versus  phase 
(averaging  over  100  periods,  25  phase  sectors) 


The  accuracy  of  the  averaging  process  can  be  improved  by 
computing,  for  each  phase  sector,  the  average  phase  as  well  as 
the  average  signal.  This  procedure  is  particularly  useful  in  the 
cases  when  the  signal  has  a  complex  dependence  on  the  phase 
difference,  as  in  the  examples  shown  below. 

Finally,  since  the  interaction  is  energy  dependent  and  the 
gain  factor  of  the  amplifier  fluctuates,  the  data  are  grouped 
according  to  the  input  energy  as  recorded  by  the  detector  . 

:V.4.  Nonlinear  Transmission 

With  the  heat  pipe  cold,  the  detectors  in  Do  and  D2 
record  the  interferometric  autocorrelation  shown  in  Fig.  4.  The 
data  acquisition  microcomputer  stores  also  the  average  of  this 
function  over  each  period,  which  is  the  intensity  autocorrelation 
shown  in  Fig.  7a.  Depending  on  the  tuning  of  the  laser  and  its 
peak  power,  a  variety  of  deformation  of  the  transmitted  auto¬ 
correlations  is  observed.  For  instance,  in  the  recording  of  Fig. 
7b,  taken  for  a  pressure  of  10  torr  of  lithium  vapor,  it  appears 
as  if  the  top  half  of  the  intensity  autocorrelation  has  been 
"chopped  off".  A  similar  change  is  observed  for  the  interfero¬ 
metric  autocorrelation. 
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Figure  7  Fringe  averaged  recording  on  D3  (a)  without  lithium  (top 
figure)  and  (b)  with  10  torr  of  lithium  (lower  figure). 


The  tuning  characteristics  of  the  nonlinear  transmission  are 
illustrated  by  the  computer  simulations  of  Fig.  8,  which  show 
successive  interferometric  (top)  and  fringe  averaged  (bottom) 
recordings  of  the  second  harmonic  of  the  transmitted  signal 
versus  delay,  for  various  detunings.  The  pulse  "area"  (the  time 
integral  of  the  generalized  Rabi  frequency,  which  is  proportional 
to  the  pulse  energy  density)  is  equal  to  1.  The  calculations  are 
carried  out  in  the  "thin  sample"  approximation  (negligible 
propagation  effects).  Far  below  resonance  (right  side  of  Fig.  8, 
'=  12  1011  s’1,  the  pulse  doesn't  interact  with  the  line,  and  the 
autocorrelations  are  unchanged  from  the  ones  without  lithium.  At 
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a  detuning  of  9  10 x x  s " x ,  which  corresponds  to  6  inverse  pulse- 
widths  below  resonance,  the  top  of  the  autocorrelation  is 
truncated.  This  is  because  near  zero  delay,  the  stark  shift 
corresponding  to  the  coherent  superposition  of  the  two  pulses  is 
sufficiently  large  to  bring  the  resulting  peak  electric  field 
into  resonance.  This  corresponds  also  to  the  situation  shown  in 
the  experimental  recording  of  Fig.  7. 
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Figure  8:  successive  Interferometric  (top)  and  Fringe  Averaged 

traces  of  the  second  harmonic  of  the  transmitted  signal 
versus  delay,  for  various  detunings  (calculated). 


With  increasing  laser  frequency  (decreasing  detuning),  the  "hole" 
burnt  by  the  resonance  line  in  the  autocorrelations  moves  to 
larger  delays.  In  the  last  case,  the  coherence  region  extends 
very  far  in  the  wings  of  the  interferometric  autocorrelation,  far 
beyond  the  limits  of  pulse  overlap.  Near  zero  delay,  the  second 
harmonic  of  the  transmitted  signal  can  be  larger  than  that  of  the 
incident  pulses,  because  of  pulse  reshaping  (shortening). 

The  third  harmonic  intensity  varies  also  with  delay, 
frequency  and  phase,  The  maxima  (for  each  27rphase  interval)  of 
uhe  third  harmonic  varies  wi t h  delay  as  shown  by  the  dashed 
curves  in  Fig.  8.  The  envelope  of  the  third  harmonic  versus  delay 
far  below  resonance  (Ivj  =  12  lO-'-1  s  ~  ■*• )  is  a  third  order 

interferometric  autocorrelation,  with  a  peak  to  background  of  32 
to  1.  In  general,  the  shape  of  the  upper  envelope  of  the  third 
harmonic  versus  delay  follows  qualitatively  that  of  the  second 
harmonic  signal  versus  delay,  with  the  largest  harmonic 
conversion  being  generally  close  to  zero  delay.  The  situation  is 
different  in  the  case  of  propagation  through  thick  samples.  For 
an  interpulse  delay  between  5.5  and  8  psec,  at  3  1011  s“l  below 
resonance  (region  indicated  by  vertical  lines  in  Fig.  8),  the 
anomalous  transmission  of  the  "90°  phase  shifted"  pulse  sequence 
exceeds  by  one  order  of  magnitude  the  "in  phase"  sequence, 
because  of  a  near  lossless  propagation  of  the  second  pulse 
Phased  matched  third  harmonic  generation  should  reach  peak 
efficiencies  (total  energy  conversion  efficiency  over  the  entire 
beam  cross  section)  of  the  order  of  1%  in  these  conditions. 


IV-5.  Data  on  Nonlinear  Transmission  and  Harmonic  Generation 

The  data  of  nonlinear  transmission  and  harmonic  generation 
versus  delay,  taken  for  various  detuning  conditions,  show  the 
characteristic  patterns  of  Fig.  8.  As  in  the  theoretical 
simulations,  as  the  detuning  (laser  below  resonance)  is 
increased,  the  "dip"  in  the  upper  envelope  of  the  inter ferogram 
of  the  transmitted  signal  (in  second  harmonic)  moves  towards  the 
center  of  the  autocorrelation  (for  that  detuning,  only  the  sum  of 
the  two  pulses  at  zero  delay  has  a  sufficient  Stark  shift  to 
"reach"  the  resonance).  A  typical  recording  is  shown  in  Fig.  9. 
The  dashed  line  ( t r i angles )  is  a  recording  of  the  interferometric 
autocorrelation  of  the  pulses.  The  dotted  (squares)  and  solid 
(circles)  lines  in  Fig.  9  are  the  envelope  of  the  recording  of 
the  transmitted  signal  (in  second  harmonic)  and  of  the  third 


o 


Figure  9 

Envelope  of  the  incident  (dashed  lines) 
and  transmitted  signal  (dotted  lines) 
and  of  the  third  harmonic  (solid  lines). 

harmonic,  respectively.  The  "dip"  at  a  delay  of 
corresponds  to  the  delay  at  which  the  Stark  shift  due 
coherent  superposition  of  the  fields  equals  the  detuning. 
Another  recording  taken  closer  to  resonance  is  shown  in  F 


Figure  10 
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Consistently  with  the  numerical  calculations,  the  transmitted 
signal  near  zero  delay  appears  to  be  larger  than  the  incident 
signal.  This  apparent  amplification  is  rnerily  the  manifestation 
of  pulse  reshaping  (compression)  near  zero  delay.  The  second 
harmonic  of  the  (shorter)  transmitted  intense  pulses  near  zero 
delay  is  larger  than  the  second  harmonic  energy  contained  in  the 
incident  pulses. 


IV. 6.  Phase  Dependence 


Figure  11  shows,  for  a  particular  inter- pulse  delay  of  2 
psec,  the  fundamental  (measured  in  second  harmonic)  before 
(dashed  line)  and  after  (dotted  line)  the  lithium  heap  pipe,  as 
•.■ell  as  the  third  harmonic  (solid  line)  as  a  function  of  the 
relative  phase  between  two  pulses.  The  phase  interval  has  been 
divided  in  16  sectors.  The  recordings  of  Fig.  11  represent  the 
average  of  1  GO  successive  periods.  The  data  points  at  a  delay  of 
2  psec  in  Fig.  9  represent  the  maxima  and  minima  of  the  curves  of 
iig.  11.  For  this  relatively  small  delay  (compared  to  the  pulse 
duration),  the  phase  dependence  is  dominated  by  the  constructive- 
destructive  interferences  of  the  overlapping  pulses  at  the 
fundamental  wavelength.  For  large  delays,  the  second  pulse  is 
interfering  constructively  and  destructively  with  the  multiphoton 
excitation  created  by  the  first  one,  in  addition  to  interfering 
partially  with  itself  (in  the  wings  of  the  pulse).  At  phase 
shifts  of  90°  and  270  °,  the  second  pulse  is  amplified  by  two- 
photon  stimulated  emission,  resulting  in  an  enhanced 
transmission.  For  delays  exceeding  the  pulse  duration,  the 
transmitted  signal  versus  phase  has  a  periodicity  of  tt. 


Figure  11 
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? i g u r a  12:  Transmitted  signal  and  third  harmonic 
versus  phase  for  a  delay  of  2.5  psec. 

At  intermediate  delays  such  as  in  the  case  of  Fig.  12,  the 
transmitted  signal  (dashed  line)  is  still  dominated  by  the 
constructive-destructive  interferences  in  the  region  of  pulse 
overlap.  However,  the  third  harmonic  shows  the  predicted 
enhancement  at  dephasings  of  90°  and  270  °.  There  is  another 
interesting  feature  in  Fig.  12,  that  also  appears  in  all 
measurements  corresponding  to  large  interpulse  delays.  The 
periodic  third  harmonic  signal  is  dephased  with  respect  to  the 
transmitted  signal.  This  dephasing  angle  is  -77/6,  and  corresponds 
to  the  phase  angle  of  the  complex  third  order  susceptibility.  We 
have  found  indeed  in  previous  theoretical  work  1  the  third 


harmonic  generation  per  unit  "characteristic  distance' 
proportional  to  the  quantity:  r  ..  .  „  a 

The  ir. aginary  part  of  thi 
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nonlinear  coefficient  -  hence  the 
,o  -  is  due  to  the  contribution  of  the 
continuum.  The  measurement  of  the  phase  dependence  of  the  third 
harmonic  offers  thus  a  unique  way  to  measure  the  contribution  (to 
the  “hird  order  susceptibility)  of  the  continuum  relative  to  that 
of  t h o  resonance  ,  or  more  generally,  the  phase  angle  of  thf 
nonlinear  suscept ibi 1 ity. 
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IV. 7.  phase  Matching. 


Theoretical  calculations  ^  predicted  a  ratio  of  2  to  1  for  a 
phase  matched  lithium-magnesium  mixture  in  conditions  of  coherent 
propagation  (a  ratio  completely  different  from  that  corresponding 
to  "incoherent"  or  "small  signal"  phase  matching).  The 
transmitted  signal  versus  delay  and  the  envelope  of  the  third 
harmonic  versus  delay  in  a  2/1  mixture  of  magnesium/lithium  are 
shown  in  rig.  13.  A  maximum  conversion  efficiency  of  6%  was 
predicted  for  an  infinite  plane  wave,  propagating  through  at 
least  four  t wo-photon  characteristic  distances  (one  characte¬ 
ristic  distance  is  1.4  cm  at  the  operating  pressure  of  10  torr). 
The  maximum  achievable  conversion  efficiency  over  the  entire 
cross-section  of  a  finite  beam  is  of  course  smaller  by  approxima- 
taly  one  order  of  magnitude.  A  conversion  efficiency  between  0.5 
end  11  is  measured.  Even  though  these  numbers  appear  very 
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Figure  13 

Third  Harmonic  and  transmitted  signals  versus  delay. 
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close  to  the  theoretical  limit,  it  should  be  noted  that  the 
measured  efficiency  is  still  limited  by  the  power  available  from 
the  source.  Indeed,  the  confocal  parameter  of  the  beam  focused  in 
the  heat  pipe  is  only  10  cm,  hence  too  short  for  the  requirement 
of  uniform  beam  cross  section  over  that  distance.  The  energy 
density  in  vacuum  varied  in  this  distance  from  50  mJ/cm^  to  100 
mJ/cm  (while  the  thick  sample  calculation  call  for  a  constant 

y 

energy  density  of  100  mJ/cm^  over  that  distance).  The  similarity 
between  fig.  13  and  the  "thin  sample"  simulations  of  Fig.  3  is 
due  to  the  too  short  confocal  parameter  of  the  focusing  in  the 
sample.  An  additional  amplification  stage  would  be  required  to 
reach  conversions  over  1%,  peaking  at  delays  larger  than  2  psec. 

I V . 3 .  Conclusions 

Enhancement  of  third  harmonic  generation  through  coherent 
interactions  in  conditions  of  two  photon  resonances  has  been 
demonstrated  for  the  first  time  experimentally.  The  agreement 
between  theoretical  prediction  and  the  experimental  results 
indicates  that  self  focusing  and  defocusing  effects  are  not 
important  enough  to  affect  the  conversion  efficiency. 

This  series  of  experiments  demonstrates  that  it  is  possible 
to  construct  a  picosecond  dye  amplifier  with  a  gain  over  106 
without  introducing  spatial  (wavefront)  or  temporal  (chirp)  phase 
distortion.  Phase  distortion  due  to  convection  currents  in  the 
1000°C  heat  pipe  were  also  found  to  be  negligible. 

Recording  of  the  third  harmonic  and  transmitted  signal 
versus  phase  enabled  us  to  measure,  for  the  first  time,  the  phase 
angle  of  the  third  order  susceptibility  responsible  for  third 
harmonic  generation.  The  corresponding  imaginary  part  of  the 
nonlinear  susceptibility  was  seen  to  match  exactly  the  calculated 
contribution  of  the  continuum. 

The  purpose  of  the  experiment  in  lithium  was  to  demonstrate 
the  points  cited  above ,  which  can  be  transposed  to  any  scheme  of 
up-conversion  in  presence  of  a  two-photon  resonance.  The  case  of 
a  higher  order  (four  photon)  resonance  is  considered  in  the  next 
section . 


V  COHERENT  INTERACTIONS  IN  MERCURY  VAPOR 
V. 1 .  Introduction 

From  the  theoretical  point  of  view,  the  case  of  a  two  photon 
resonance  is  a  very  particular  case  of  interaction  where  all 
terms  are  second  order  in  the  field.  For  a  negligible  Stark  shift 
and  photoionization,  there  is  only  one  parameter  (apart  from  a 
scaling  factor)  distinguishing  one  atom  from  another:  the  ratio 
of  the  probability  of  third  harmonic  generation  to  that  of  two- 
photon  stimulated  emission  (£i2//'*i?  in  the  notations  of  ref.  1). 
There  is  no  optimal  pulse  duration,  within  the  constraint  that 
the  pulse  duration  be  shorter  than  the  elastic  collision  time. 

The  case  of  higher  order  resonances  (N)  is  much  more 
complex:  the  equations  will  contain  terms  of  all  order  up  to  the 

The  case  of  four  photon  resonant  excitation  of  the  6s  -  7s 
transition  of  mercury  vapor  is  particularly  interesting  but 
complex,  because  of  a  near  three  photon  resonant  enhancement 
through  the  6p  level.  A  theoretical  investigation  of  coherent 
propagation  and  harmonic  generation  in  mercury  vapor  has  been 
completed,  using  the  same  formalism  as  for  the  study  of  harmonic 
generation  in  lithium  vapor  Preparations  for  an  experiment  to 
complement  the  theory  have  been  made  (laser  source,  vapor  cell). 

V . 2 .  Theory 

As  in  the  case  of  a  two  photon  resonance,  a  generalized  Rabi 
frequency  can  be  defined,  which  is  here  proportional  to  the 
fourth  power  of  the  electric  field  of  the  light  at  the 
fundamental  frequency.  However,  the  analogy  stops  here,  because 
the  presence  of  terms  of  lower  order  makes  it  impossible  to 
express  all  electric  field  related  quantities  in  units  of  the 
Rabi  frequency. 

The  equations  describing  coherent  four-photon  resonant 
propagation  and  third  and  fifth  harmonic  generation  in  mercury 
vapor  are  detailed  in  Appendix  B.  The  dipole  matrix  elements 
needed  to  evaluate  all  the  coefficients  of  these  equations  were 
calculated  from  the  absorption  coefficients  corresponding  to  the 
main  transitions.  Some  major  conclusions  from  the  theory  are 
listed  below. 


For  large  delays,  the  phase  dependence  of  the  transmitted 
signal,  the  third  and  fifth  harmonic  have  a  periodicity 
of  7T/2. 

For  picosecond  pulses  at  560  nm,  an  energy  density  of  the 
order  of  1  j/cm2  is  required  to  significantly  populate 
the  upper  level  of  the  resonant  transition. 

For  femtosecond  pulses  at  620  nm,  the  energy  density 
required  to  pump  significantly  the  7s  level  is  of  the 
order  of  100  mJ/cm2. 

The  proximity  of  the  6p  level  at  3/4  of  the  6s  -  6d 
transition  results  in  an  anomalously  large  Stark  shift 
and  an  enhancement  of  third  harmonic  generation! 

A  large  photo  ionization  cross-section  from  the  6d  level 
would  prevent  efficient  conversion  to  higher  harmonics  in 
a  two  pulse  sequence. 

The  theory  leaves  several  basic  questions  unanswered,  that  will 
have  to  be  resolved  by  experiments.  First,  the  Stark  shift  of  the 
6s  -  6d  transition  (4  photons  at  560  nm)  is  so  large,  that  high 
energy  femtosecond  pulses  at  620  nm  could  be  simultaneously 
resonant  with  the  6s  7s  and  with  the  6s  6d  transitions!  While  it 
is  likely  that  the  theoretical  model  chosen  breaks  down  in  that 
situation,  there  might  still  be  a  large  enhancement  of  the  non¬ 
linear  susceptibilities. 

Another  basic  unanswered  question  pertains  to  the  large  photo¬ 
ionization  from  the  excited  level.  Has  the  photoelectron  lost 
complete  memory  of  its  position  in  the  Hg  atom,  even  10  fsec 
after  its  ejection?  This  question  can  be  tested  by  counting  the 
number  of  ions  created  by  a  pair  of  closely  spaced  pulses,  as  a 
function  of  their  relative  phase. 
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V . 3 .  Source  Developments 


For  a  given  pulse  energy,  the  shorter  pulses  of  higher  peak 
intensity  will  be  the  most  efficient  in  exciting  a  multiphoton 
resonance  of  order  higher  than  2.  We  have  therefore  started  the 
development  of  a  femtosecond  laser  oscillator  to  replace  the  6 
psec  synchronously  mode-locked  laser.  That  laser  oscillator 
combines  the  features  of  the  passively  mode-locked  dye  laser  ' 
with  the  tunability  of  the  synchronously  mode-locked  laser. 
Passive  mode -locking  is  achieved  in  an  "antiresonant  ring" 
located  at  one  end  of  the  cavity  10.  The  antiresonant  ring  is 
mounted  on  a  translation  stage,  which  allows  continuous  tuning  of 
the  cavity  length,  an  essential  feature  for  synchronous  mode¬ 
locking  2 2 .  In  recent  tests,  passive  mode-locking  alone  has 
provided  highly  stable  operation  with  50  fsec  pulses,  and  an 
average  output  power  of  50  mW  x  - .  Geometrical  astigmatism 
compensation  is  used  to  maximize  the  intracavity  power  density  in 
the  saturable  absorber  12  (Appendix  C).  The  pulse  shape  and  phase 
modulation  are  accurately  characterized,  using  techniques  of 
pulse  shape  and  chirp  determination  that  we  recently  developed 
The  "downchirp"  of  the  pulses  can  be  adjusted  to  yield  pulse 
compression  instead  of  pulse  broadening  through  the  normally 
dispersive  amplifier  chain.  A  hybrid  passively-actively  mode- 
locked  laser  based  on  this  design  is  presently  being  tested.  It 
is  hoped  that  the  hybrid  laser  will  provide  a  higher  output  power 
(100  mW)  and  easier  tunability  than  the  purely  passively  mode- 
locked  laser.  It  should  be  noted  however  that  the  present  tunabi¬ 
lity  range  of  the  purely  passively  mode-locked  antiresonant  dye 
laser  (605  to  630  nm)  is  quite  adequate  to  study  coherent  four 
photon  resonant  interaction  on  the  6s  -  7s  and  on  the  (Stark 
shifted)  6s  -  6d  transitions  of  mercury.  The  accurate  knowledge 
of  the  pulse  shape  and  modulation^  is  a  particularly  important 
asset  in  the  study  of  coherent  effects  in  photoionization. 

The  use  of  shorter  pulses  makes  it  possible  to  use  much 
higher  pressures  in  the  sample  cell  (since  the  requirement  that 
the  collision  rate  be  slower  than  the  inverse  pulse  width  is 
relaxed  by  two  orders  of  magnitude).  The  possibility  of  using 
higher  pressures,  smaller  confocal  parameters,  alleviates  much  of 
the  difficulties  associated  with  the  design  of  a  sample  cell.  A 
short  mercury  heat  pipe  has  been  constructed,  in  which  two 
stainless  steel  cones  serve  the  dual  purpose  of  ion  counting 
electrodes  and  return  channel  for  the  condensing  mercury. 


VI  USE  OF  MULTIPLE  PULSE  SEQUENCES 

We  have  extended  the  theory  to  the  use  of  multiple  pulse 
sequences  in  the  study  of  multi photon  ionization  of  lithium.  The 
results  of  this  work  are  contained  in  a  paper  published  in  Phys 
Rev.  1  ^  and  included  as  Appendix  D.  This  theoretical  study  shows 
how  random  pulse  sequences  should  be  used  to  study  the  influence 
of  source  coherence  on  a  multiphoton  process.  Another  basic  point 
arising  from  this  theory  is  that  regular  pulse  sequences  can 
provide  spectral  selectivity  well  beyond  the  pulse  bandwidth.  For 
instance,  by  tuning  the  average  frequency  of  a  train  of  10 
pulses,  the  number  of  ions  produced  in  lithium  vapor  (condition 
of  two-photon  resonance)  versus  laser  frequency  displays  a  series 
of  peaks  reminiscent  of  the  Fourier  spectrum  of  the  excitation. 
The  width  of  these  individual  pulses  can  be  used  as  a  measure  of 
the  photoionization  rate  from  the  excited  (4s)  level. 

The  basic  conclusion  from  this  theory  is  that  the  use  of 
pulse  sequences  offers  a  solution  to  the  conflicting  requirements 
of  spectral  selectivity  and  high  peak  intensities  needed  to  study 
multiphoton  processes.  This  conclusion  applies  to  molecules  as 
well  as  to  atoms.  The  study  has  been  recently  extended  to  multi- 
photon  vibrational  spectroscopy.  Pulse  sequences  have  been 
designed  to  excite  to  a  high  energy  level  an  anharmonic  ladder  of 
levels.  The  relative  phase  between  successive  pulses  determines 
the  direction  of  flow  of  energy  within  the  molecule:  whether 
mere  will  be  absorption  of  emission  of  a  photon  within  the 
ladder  of  levels.  Calculations  of  multiphoton  excitation  of  an 
anharmonic  ladder  of  vibrational-rotational  levels  have  been 
made.  The  model  chosen  was  that  of  the  molecule  CH^F.  The  calcu- 
lations  show  -L“'  -L-)  that  simple  four  or  five  pulse  sequences  can 
be  designed  to  bring  most  of  the  population  in  the  level  v  =  5. 
These  particular  pulse  sequences  result  in  a  cooling  of  the 
rotational  temperature  in  the  excited  state,  pulse  sequences 
offer  thus  the  possibility  of  achieving  vibration ally  as  well  as 
rotationally  selective  multiphoton  vibrational  excitation  of 
molecu  les  . 
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t.  Introduction 

We  report  -ho  ckveloptiu  r.t  <>t  a  pi</o>tJc»:id  pulse 
dye-laser  oscillator -amplifier  system  designed  for  ex¬ 
periments  i>t  resonant  coherent  propagation.'  Such 
experiments  put  stringent  requirements  on  the  prop¬ 
erties  of' the  laser  system,  most  of  which  will  he  reviewed 
briefly  in  Sec.  II.  We  have  been  able  to  meet  these  re¬ 
quirements  using  a  commercially  available  synchro¬ 
nously  pumped  dye  laser  in  combination  with  a  home- 
built  amplifier.  As  we  will  explain  in  Sec.  III.  this  am¬ 
plifier  consists  of  three  stages  being  pumped  by  a  fre¬ 
quency-doubled  N’dtVAd  laser.  I  he  first  two  stages 
are  pumped  transversely,  while  the  last  stage  is  pumped 
longitudinally.  We  suceeded  in  achieving  a  high  gain, 
while  avoiding  a  buildup  of  amplified  spontaneous 
emission  (ASK*,  without  the  use  of  spatial  filters  or 
saturable  absorbers.  Accurate  synchronization  of  the 
various  pulses  combined  with  a  careful  optimization  of 
the  pump  and  amplified  beam  geometry  ensures  that 
the  energy  of  the  pump  is  channeled  into  the  picosecond 
pulse  rather  than  into  ASK.  Compared  with  purely 
longitudinally  pumped  systems,  our  transverse-longi¬ 
tudinal  pumping  geometry  combines  the  advantage  of 
a  higher  output  energy  with  the  ability  to  control  the 
spatial  beam  profile.  In  Sec.  IV,  we  will  present  inter¬ 
ferometric  autocorrelation  measurement performed 
for  the  first  time  on  such  a  picosecond  pulse  dve-laser 
oscillator  amplifier  system.  These  measurements,  for 
which  a  new  kind  of  data  acquisition  system  needed  to 
he  developed,  show  that  our  amplifier  does  not  signiti- 


The  authors  .i r»*  with  North  Trvis  State  (  imcrsitv.  (Yntor  for 
<1  ijti.intum  Kirctronics.  l*h\Mt  -  Prp.irtment.  Ikntnn.  Texas 

7b. !« 

Received  |f>  November  l!WH. 

♦  H  *»  ’  U.'.J.  I  yjloH-rd.slt'J  Mp  0. 

■  1  n  I  Opt  ical  Society  ot  America. 


cantly  change  the  temporal  properties  of  the  dye-laser 
pulses.  Preliminary  measurements  of  coherent  two- 
photon  absorption  in  lithium  demonstrate  that  this 
system  meets  the  stringent  coherence  requirements  of 
a  multiphoton  coherent  propagation  experiment. 

II.  General  Design  Considerations  for  a  Picosecond 
Pulse  Dye  Amplifier 

Even  though  we  designed  our  amplifier  for  achieving 
a  particular  goal  (saturation  of  the  2s -4s  transition  of 
lithium  at  571.2  nmi,  the  main  specifications  are  de¬ 
sirable  requirements  for  a  large  number  of  experiments: 
(1)  high  output  energy;  (2*  minimum  ASK;  C>)  good 
beam  quality  for  interferometric  experiments,  i.e.. 
uniform  beam  profile  and  Hat  distortion-free  wave 
fronts;  (4 >  no  temporal  phase-modulation  efferts;  (.*>) 
both  short-  (pulse-to-pulsel  and  long-term  stability;  tdl 
availability  of  a  nv  pilot  beam  propagating  through  the 
system  for  the  ease  of  alignment  of  the  target  experi¬ 
ment.  Each  requirement  dictated  one  or  more  basic 
design  features  of  the  amplifier  to  he  built  as  we  will 
explain  now.  (See.  e.g..  Refs.  4  and  5  and  additional 
references  therein  fora  theoretical  analysis.) 

Obviously,  a  high  output  energy  calls  for  an  optimum 
use  of  the  pump  radiation,  which  in  turn  requires  the 
minimum  beam  diameter  in  each  amplifier  stage  com¬ 
patible  with  a  maximum  unsaturated  gain.  Maximizing 
the  gain  for  a  given  input  also  requires  minimization  of 
ASK.  This  can  be  achieved  through  1 1 )  accurate  timing 
of  the  pump  pulse  relative  to  the  picosecond  pulse  to  be 
amplified  in  each  stage,  thereby  minimizing  the  overlap 
between  the  gain  puKe  in  one  cell  and  the  ASK  trom  the 
adjacent  cell;  (21  accurate  geometrical  matching  of  the 
[lump  beam  and  the  picosecond  dve-laser  beam  to 
maximize  uniform  pump  deplot  i<  n  by  t  he  latter  in  each 
i  ell;  t:t)  appropriate  choice  ol  dyes.  The  choice  ot  dyes 
is  lor!  her  limited  by  mir  requirement  that  the  system 
should  be  transparent  for  the  unamplilied  dye-laser 
pulses,  which  serve  as  a  cw  pilot  beam.  This  last  out- 
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dition  is  particularly  important  in  applications  which 
require  a  critical  alignment  of  the  amplified  beam.  For 
our  applications,  the  dye-laser  beam  is  sent  first  along 
the  axis  of  an  interferometric  delay  line  and  subse¬ 
quently  through  a  long  but  narrow  heat  pipe,  an  oper¬ 
ation  that  requires  a  cw  beam.  The  requirement  of  an 
aberration-free  temporal  and  spatial  beam  profile  im¬ 
poses  the  following  design  features:  ( 1 )  short  dve  cells 
to  minimize  pulse-broadening  and  phase-modulation 
effect.-*  due  to  the  solvi  nt1’-':  C_)  the  absence  of  Filters  and 
saturable  absorbers  to  avoid  wave  front  distortion  in 
subsequent  amplifiers;  to)  a  choice  of  solvents  to  obtain 
maximum  thermal  conductivity;  (4)  careful  lateral  ad¬ 
justment  of  the  pump  beam  in  the  longitudinally 
pumped  cell  to  properly  shape  the  spatial  beam  profile 
ot  the  amplified  dye-laser  pulses.  Finally,  stability 
requirements  put  stringent  conditions  on  the  pump 
source  for  the  amplifier:  tl)  short-term  stability  de¬ 
pend'  on  the  synchronization  between  dye  laser  and 
amplifier  pump;  1 2)  long-term  stability  merely  imposes 
a  constant  thermal  lens  effect  for  the  pump  laser. 

III.  System  Descriptions 

The  oscillator  is  a  Spectra-Physics  dye  laser  syn- 
chronousiy  pumped  by  a  Spectra-Physics  Ar-ion  laser 
actively  mode  locked  at  82  MHz.  For  the  gain  medium 
we  use  rhodamine  6G  in  ethylene  glycol  (concentration 
2  X  lo~  *  NI).  Although  rhodamine  560  would  yield  a 
slightly  higher  gain  at  571  nm.  we  prefer  instead  to  use 
rhodamine  6(1  because  of  its  better  chemical  stability. 
(The  rhodamine  560  solution  has  to  be  replaced  every 
week.)  Wavelength  tuning  ~571  nm  is  achieved  with 
a  three-plate  hirefringent  filter.  With  a  50%  output 
coupler,  this  sytem  generates  a  continuous  train  of 
5.6-nJ  pulses  with  a  duration  of  6.2-psec  FWHM. 

Our  three-stage  dye  amplifier  sketched  in  Fig.  I  is 
pumped  by  a  frequency-doubled  (^-switched  Nd:YAG 
oscillator-amplifier  (Quanta  Ray  DCR-1).  This  system 
delivers  280-mJ  pulses  (560  nm)  of  8-nsec  FWHM  at  a 
repetition  rate  of  10  Hz.  The  polarization  of  this  beam 
is  orthogonal  to  the  polarization  of  the  dve-laser  beam. 
Therefore,  we  put  the  longitudinally  pumped  stage 
between  crossed  polarizers  to  couple  out  both  the  am¬ 
plified  dye-laser  beam  and  the  remaining  part  of  the 
pump  beam.  A  1.5-rn  optical  delay  separates  the  dye- 
laser  oscillator  from  the  first  amplifying  cell  to  prevent 
oscillation  of  the  latter  with  the  output  mirror  of  the 
former. 

A  Beam  Geometry 

To  prevent  ASK  from  building  up  between  stages,  a 
direct  vision  prism  is  inserted  between  the  first  and 
second  stages.  This  is  the  sole  optical  isolator  used  in 
our  amplifier,  in  contrast  to  other  systems  which  need, 
in  addition.  >aturab!e  absorbers  and  pinholes.4 'y 
1'  iriicnlarly  important  in  this  respect  i'  the  expanding 
beam  geomet ry  throughout  the  system  not  only  to  in- 
erca.'f  t  he  effect  iveness  of  the  prism  in  discriminating 
against  ASK  hut  also  to  maximize  the  overall  gain  of  the 
amplifier.  Starting  with  a  diameter  of  <1  mm  in  the 
lif't  cell,  the  beam  diameter  is  increased  three  times  in 


Fig  1.  Schematic  representation  of  the  three-stage  amplifier  (DBS 
=  tiichroic  beam  splitter,  POL  =  polarizing  beam  splitter,  P  =  direct 
vision  prism). 


the  second  cell  and  further  expanded  in  the  last  stage. 
The  latter  is  inserted  between  two  lenses  of  a  collimating 
beam  expander  with  an  aperture  ratio  of  5.  This  design 
configuration  of  the  last  stage  serves  a  dual  purpose. 
The  expansion  of  the  dve-laser  beam  is  required  to  keep 
the  energy  density  at  571  nm  below  the  saturation  en¬ 
ergy  density.  On  the  other  hand,  the  converging  ge¬ 
ometry  of  the  pump  beam  is  needed  to  maintain  a  suf¬ 
ficient  level  of  excitation  despite  absorption  in  the  dye. 
One  should  bear  in  mind  that  the  performance  limits 
of  the  dye-amplifier  are  defined  in  the  last  stage.  In¬ 
deed,  the  saturation  energy  density  over  the  output 
cross  section  of  the  last  stage  is  the  maximum  energy  per 
pulse  that  the  system  can  deliver  (e.g.,  1  mJ/cm2  for 
rhodamine  6G).KI  Because  of  its  large  cross  section,  this 
last  stage  requires  most  of  the  energy  (77%)  of  the  pump 
laser,  while  yielding  only  a  low  gain  (factor  8).  Of  the 
remaining  pump  energy.  7%  is  used  for  the  first  stage 
(gain:  10!)  and  16%  for  the  second  stage  (gain:  50). 

B.  Choice  of  Dyes 

For  the  transversely  pumped  stages,  we  use  1-cm  long 
flow  through  cells  at  Brewster  incidence.  The  best  dye 
candidate  for  a  pump  wavelength  near  550  nm  is  rho¬ 
damine  6G  (highest  gain  at  571  nm,  longest  lifetime). 
However,  the  absorption  hand  of  rhodamine  6G  extends 
well  beyond  571  nm.  Therefore,  rhodamine  560  has  to 
be  used  in  the  First  two  stages  (where  the  concentration 
is  necessarily  high)  to  obtain  a  reasonably  large  trans¬ 
mission  of  the  un  unplified  dye-laser  pulses  through  the 
system.  It  was  determined  empirically  that  the  fol¬ 
lowing  dye  combinations  in  the  first  and  second  cells 
yield  the  best  results  (maximum  pump  depletion, 
minimum  ASK.  and  acceptable  loss  at  571  nm):  for  the 
first  stage,  a  mixture  of  rhodamine  560  and  rhodamine 
6G  in  methanol  (respectively,  2  X  10- 1  M  and  10-4  Mt; 
for  the  second  stage,  rhodamine  560  in  methanol  ( 10“  ' 
M).  Finally,  for  the  last  stage,  which  is  5  cm  long,  we 
use  rhodamine  6G  to  fullly  maximize  the  gain.  The 
solvent  in  this  stage  is  a  <'.5%  aqueous  solution  ot  Am- 
monix  L.O.  to  minimize  thermal  effects.  For  the  small 
dye  concentration  in  this  final  stage  ( 1 .5  X  lt)“;’  M),  the 
absorption  of  the  unamplified  pulses  is  negligible. 
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C.  Synchronization 

In  Sir.  ]1  we  emphasized  the  need  !<>r  accurate  and 
jit tt  r-l rt-e  >vnchronization  in  each  Maui*  between  the 
pump  pulse  and  the  dye-laser  pulse  to  be  amplified.  To 
reneh  this  „oai,  v.e  use  the  mode-locker  drive  of  the 
Ar-ion  laser  as  an  absolute  dock.  Through  direct  dig¬ 
ital  division,  this  clock  signal  is  staled  down  to  ^10  Hz. 
Next,  this  signal  is  used  to  interrogate  the  original  dock 
•  agnal  in  a  logical  \M>  gate  thus  providing  a  in-Hz.  ref¬ 
erence  signal  which  is  synchronized  with  the  master 
clock  (to  within  better  than  O.n  nsect.  A  detailed  de¬ 
scription  of  the  electronic  circuit  lias  been  published 
elsewhere.’1  The  lu-Hz  reference  in  turn  is  used  to 
directly  trigger  the  Q-switch  and.  after  an  appropriate 
delay,  also  the  rtashlamps  of  the  N'd:YA(»  oscillator- 
amplifier.  The  resulting  jitter  in  the  synchronization 
between  the  pump  pulse  and  the  dye-laser  pulse  to  he 
amplified  is  of  the  order  of  1  nsec. 

To  further  minimize  energy  fluctuations  of  the  am¬ 
plified  pulse  due  to  this  small  remaining  jitter,  we  syn¬ 
chronize  the  first  stage  at  the  leading  edge  of  the  pump 
pulse,  while  the  second  stage  is  synchronized  at  the 
trailing  edge.  This  procedure  further  helps  in  reducing 
undesirable  coupling  in  ASH  between  the  initial 
stages. 

D  Performance  of  the  System 

The  output  energy  of  our  amplifier  is  ~l.a  m-J,  pulse. 
Pulse-to-pulse  stability  is  better  than  20r<.  This  is  a 
good  result  for  a  chain  of  amplifiers,  none  of  which  is 
saturated.  I  sing  fresh  dyes,  the  long-term  stability  is 
also  good  i  Wr  over  a  period  of  10  h).  For  an  optimum 
-y uehronization.  the  ratio  of  ASK  energy  to  pulse  energy 
after  the  la>t  stage  is  only  a  few  percent. 

The  spatial  intensity  distribution  of  the  amplified 
beam  has  also  been  examined  using  an  optical  multi¬ 
channel  analyzer  fOMAf.  Because  the  firsi  two  stages 
are  pumped  transversely,  one  does  not  expect  a  uniform 
gain  profile  through  the  whole  chain.  Therefore,  one 
may  tear  that  the  amplifier  design  may  m>t  provide  a 
useful  and  smooth  pulse  spatial  profile.  Our  mea¬ 
surements  instead  have  shown  that,  hv  carefully  po-i- 
tioning  the  pump  beams  relative  to  the  amplified  beam 
in  the  last  two  stages,  we  can  achieve  complete  control 
of  f  he  prolile  ol  the  output  beam.  An  example  is  shown 
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in  Fig.  2.  With  the  OMA  in  the  tree-running  mode,  no 
shot-to-shot  variations  of  this  profile  were  detected. 

The  quality  of  the  amplifier  output  was  also  exam¬ 
ined  with  a  Mach-Zehnder  interferometer.  Complete 
extinction  of  the  beam  over  its  entire  cross  section  could 
he  achieved  with  the  interferometer  balanced.  This 
demonstrates  the  absence  of  wave  front  distortion  in  the 
amplifier  chain. 

IV.  Temporal  Properties  of  the  Amplified  Pulses 

Kquallv  important  for  coherent  propagation  experi¬ 
ments  are  the  temporal  properties  of  the  amplified 
pulses.  We  infer  these  properties  from  second-order 
intensity  and  fringe-resolved  autocorrelation  mea¬ 
surements  using  a  Mach-Zehnder  interferometer.- 
The  main  difficulty  involved  in  these  measurements 
stems  from  the  fact  that  the  delay  arm  of  the  interfer¬ 
ometer  has  to  he  moved  uniformly  with  absolute  sub- 
wavelength  accuracy  over  a  sufficiently  long  distance. 
Indeed,  the  requirement  that  the  motion  of  the  delay 
arm  be  continuous  without  the  discontinuities  that 
would  result  from  the  static  friction  coefficient  (being 
larger  than  the  dynamic  one)  implies  a  minimum 
scanning  speed.  On  the  other  hand,  however,  an  upper 
limit  for  this  scanning  speed  is  set  by  the  repetition  rate 
of  the  amplifier  ( 10  Hz)  and  the  requirement  that  there 
be  at  least  twenty  to  thirty  data  points  per  period  to 
resolve  individual  fringes.  Obviously,  both  require¬ 
ments  are  conflicting.  Moreover,  the  remaining 
pulse-to-pulse  instabilities  will  cause  large  fluctuations 
<>n  top  of  second -order  fringes,  thereby  washing  out  the 
fringe  contrast.  Averaging  out  t he  data  over  a  large 
number  of  laser  shots  will  only  further  reduce  this 
contrast.  1\>  solve  these  problems,  we  had  to  develop 
an  adequate  data  acquisition  system. 

A.  Data  Acquisition  System 

The  hardware  of  mir  data  acquisition  system  is 
schematically  represented  in  Fig.  T  About  l'r  of  the 
amplifier  output  is  used  as  a  reference  beam,  the  re¬ 
maining  part  is  -cut  through  the  interferometer.  Both 
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the  rtiereme  beam  and  one  of  the  output  beams  of  the 
interferometer  are  frequency  doubled  and  digitized 
after  detection  bv  photmonductive  PIN  diodes  (I’DT 
l(H)L).  The  other  output  beam  of  the  interferometer 
is  used  for  two-photon  coherent  propagation  experi¬ 
ments  m  lithium  vapor.  The  electronic  discriminator 
in  the  reference  channel  rejects  pulses  with  either  too 
low  or  too  high  US'  energy.  This  is  an  essential  feature 
for  our  coherent  two-photon  propagation  experiments 
in  lithium  vapor  which  depend  critically  on  the  pulse 
energy  (through  induced  Stark  shifts). 

The  aforementioned  problem  of  smoothing  the  in¬ 
terference  fringes  is  solved  in  the  software  of  our  data 
acquisition  system.  First,  errors  due  to  amplitude 
fluctuations  of  the  pulses  are  reduced  to  a  few  percent 
by  computing  t  he  ratio  of  the  signals  in  channels  it  and 
1 .  Next,  while  varying  the  delay  in  the  autocorrelator 
at  a  regular  speed  of  0.1  tip m/ sec,  the  resulting  fringes 
are  partitioned  over  a  large  number  of  sets  teach  set 
typically  has  100  fringes).  Per  set,  two  different  aver¬ 
aging  processes  are  performed.  On  the  one  hand,  all 
data  in  eaeh  set  are  averaged  thereby  washing  out  the 
fringes  completely.  For  a  large  number  of  sets  this  will 
vield  the  conventional  0:1  intensity  autocorrelation 
,  urve.  On  the  other  hand,  an  average  fringe  per  set  is 
constructed  by  averaging  the  data  in  a  single  periodicity 
over  the  whole  set.  An  example  is  shown  in  Fig.  4.  This 
averaging  process  permits  one  to  keep  track  of  small 
phase  shifts  over  consecutive  sets  of  fringes.  For  a  large 
number  of  sets,  it  yields  the  fringe-resolved  8:1  auto¬ 
correlation  curve.  In  this  way  one  obtains  both  auto¬ 
correlation  curves  during  a  single  experiment. 

B.  Data  Reduction 

By  fitting  the  experimental  data  for  both  the  inten¬ 
sity  and  the  fringe-resolved  autocorrelation  curves,  it 
is  possible  to  determine  accurately  t  he  temporal  prop- 
erties  of  the  amplified  pulses.  From  the  experimental 
arrangement  described  above,  we  obtain  two  recordings 
of  second  harmonic  signal  vs  delay:  the  intensity  and 
the  interferometric  autocorrelation.  The  intensity 
autocorrelation  curve  provides  a  first  estimate  of  a 
reasonable  pulse  shape  and  pulse  duration  which  is  used 
as  a  starting  point  in  attempting  to  fit  the  interfero¬ 
metric  autocorrelation  measurements.  For  a  given  trial 
function  for  the  complex  electric  field  envelope 

/•.'< :  i  =  /-.'f  i  exp1:  •in|,  <  I ) 

a  numerical  computation  is  made  of  the  interferometric 
second -order  autocorrelation 
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I  he  computations  are  organized  as  follows:  first,  at  r 
=  ' '  and  at  r  =  I, Jr  the  first-fringe  maximum  and 
minimum  are  computed:  Mibsequenth  r  is  incremented 


o2oi  xiatu 


nut  of  -'oo.jp  t  sees  secomos 

•ssuEl  »:  ztitt  of  .soup  -  a  aeea  mcbons 

Fig.  i.  K.v.im'ple  ,m  .nerage  fringe  representing  a  set  .if  I  rinses. 


by  Ax  =  1,  ;■  to  generate  the  next  pair  of  points  and  so 
on.  To  improve  the  accuracy,  an  additional  loop  gen¬ 
erates  near  eaeh  fringe-maximum  three  points  of  F(x ), 
not  more  than  two  of  them  being  located  at  the  same 
side  ol  that  maximum.  With  a  parabolic  fit  through 
these  points,  the  computer  determines  (I)  the  true  lo¬ 
cation  of  maximum  (2)  the  correct  Ax,  which  will 

he  used  to  initiate  the  next  main  loop.  Cl)  the  location 
of  the  fringe-minimum  tximn),  and  finally  (41  the  cor¬ 
responding  values  of  Ftx,.laJ,  F(rinill).  and  /(xmax). 
Once  /( r)  and  both  envelopes  of  F(x)  are  computed  in 
this  way.  the  pulse  parameters  (shape,  asymmetry, 
width,  and  chirp)  are  adjusted  until  a  good  fit  with  the 
experimental  data  can  be  reached.  Finally,  for  the  best 
fitting  pulse,  the  deviation  of  the  average  frequency 

i  (t>E'-<t)dt 

— —  (D 

J  /•-  i(»  if 

and  the  pulse  bandwidth1- 
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are  obtained. 


C.  Results:  Pulse  Characteristics  Before  and  After 
Amplification 

Figure  "><a)  shows  the  intensity  autocorrelation  curve 
for  the  amplified  pulses  (dashed  line).  The  solid  line 
represents  a  numerical  lit  using  a  Claussian  pulse  shape 
with  6.2-psec  FWHM  |i.e.,  5.27-psec  H\V(l/<-)A/|. 
Note  that  an  intensity  autocorrelation  cannot  effec¬ 
tively  distinguish  between  a  Claussian  and  a  sech  pulse 
shape. 

From  Kq.  (2),  it  is  clear  that  the  intensity  autocorre¬ 
lation  is  phase  independent  and.  therefore,  does  not 
carry  any  information  about  the  pulse  chirp.  On  the 
other  hand,  as  shown  in  Kels.  J  and  3.  the  interfero¬ 
metric  autocorrelation  is  sensitive  to  the  pulse  chirp 
'!>U ).  The  experimental  results  obtained  for  the  latter 
with  our  amplified  pulses  are  shown  in  Fig.  at  hi  (dashed 
line).  The  solid  line  in  Fig.  nth)  represents  again  the 
best  theoretical  lit:  Gaussian  pulses  (H.J-psec  FWHM) 
with  a  phase  modulation  proportional  to  the  pulse  in¬ 
tensity  (<;>  =  I.  Id  /*,'-).  For  comparison,  the  *:1  auto- 
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Kii. Kxpi  rimt  n:  il  resi;lt<  and  best  theoretical  tits  for  (at  the  in- 
u-n-ity  autoci  rrola;  ion  curve  and  i b)  the  fringe-resolved  autocorre¬ 
lation  t  urve  (dashed  lines  are  experimental  curves,  solid  lines  repre¬ 
sent  the  theoretical  tits;  for  the  dotted  line  in  (b),  see  text]. 

correlation  curve  for  the  same  Gaussian  pulses  but 
havin','  instead  a  linear  phase  modulation  [best  fit:  <i> 
=  1.0  (G'5.2 7)-|  is  also  shown  in  Fig.  5(b)  (dotted  line). 
The  data  imply  that  the  amplified  pulses  are  not  com¬ 
pletely  bandwidth  limited;  according  to  Eq.  (5),  the 
bandwidth  is  79  GHz. 

The  data  shown  in  Fig.  5(b)  are  much  more  sensitive 
to  phase-modulation  effects  than,  e.g.,  spectral  data. 
The  bandwidth  duration  product  for  the  nonlinear 
chirped  pulses  in  our  case  is  0.49  or  only  10rb  larger  than 
the  minimum  of  0. 14  for  unchirped  Gaussian  pulses. 
On  the  other  hand,  the  FWHM  of  the  upper  envelope 
of  the  interferometric  autocorrelation  curve  for  our 
pulses  is  7.2  psec  compared  with  9.5  psec  for  unchirped 
pulses;  the  lower  envelope  is  even  more  sensitive  to 
phase-modulation  effects:  its  FWHM  changes  by  al¬ 
most  5()r;  from  our  chirped  pulses  <9.5  psec;  to  the  case 
of  unchirped  pulses  1 1 1.7  psec). 

The  kind  of  chirp  we  have  detected  for  the  amplified 
pulses,  i.e.,  a  chirp  proportional  to  the  time  derivative 
of  the  pulse  intensity.  i>  typical  of  self -phase  modulation 
of  the  pulses.  An  interferometric  autocorrelation  of  the 
unamplified  dve-laser  pulses  as  shown  in  Fig.  6  clearly 
indicates  that  the  origin  of  the  chirp  is  in  the  laser  os¬ 
cillator  rathe  r  than  in  the  amplifier  chain.  Partial  chirp 
compensation  i-  taking  place  in  the  amplifier:  the  in¬ 
terferometric  autocorrelation  curve  lor  the  unamplified 
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pulses  (the  solid  line  represents  the  best  lit;  the  dashed  line  corre¬ 


sponds  to  the  best  tit  for  the  amplified  pulses). 


pulses  reveals  a  phase  modulation  proportional  to  1.84 
E-  (rather  than  1.43  E-  for  the  amplified  pulses). 

To  conclude,  our  interferometric  autocorrelation 
measurements  as  shown  in  Fig.  5(b)- have  demonstrated 
that  the  phase  distortion  introduced  by  the  amplifier 
is  small  in  the  time  as  well  as  in  the  space  domain.  A 
peak-lo-background  of  8: 1  can  only  be  obtained  with 
spatially  coherent  pulses.  Indeed,  at  zero  delay,  the 
maximum  signal  of  eight  times  the  background  requires 
that  the  two  halves  of  the  pulse  interfere  constructively 
over  their  entire  cross  section. 

V.  Conclusions 

In  summary,  we  have  developed  a  picosecond  pulse 
dye-laser  oscillator-amplifier  system  tunable  near  570 
nm.  At  a  repetition  rate  of  10  Hz  the  output  consists 
of  1.5-tnJ  pulses.  Both  the  short-  and  long-term  sta¬ 
bility  are  excellent.  Contrary  to  other  similar  systems, 
we  have  avoided  the  use  of  either  pinholes  or  saturable 
absorbers.  Therefore,  our  system  requires  fewer 
components,  which  makes  the  alignment  easier.  De¬ 
spite  the  absence  of  pinholes  and  saturable  absorbers, 
the  ratio  of  ASE  energy  to  pulse  energy  after  the  last 
stage  is  only  a  few  percent.  An  additional  advantage 
of  our  system  over  other  short  pulse  amplifiers  is  that 
the  unamplified  dye-laser  pulses  are  transmitted 
through  the  system  thus  providing  a  cw  pilot  beam. 
This  is  an  important  feature  for  applicat  ions  requiring 
a  critical  alignment.  Our  amplifier  output  has  a  uni¬ 
form  spatial  beam  profile  and  distortion-free  wave 
fronts.  The  temporal  properties  of  the  amplified  pulses 
have  also  been  investigated  using  a  computerized  data 
acquisition  system  specially  developed  for  this  purpose. 
The  pulses  are  Gaussian  (FWHM:  6.2  psec)  and 
present  a  nonlinear  chirp  (1.4  E).  The  bandwidth  is 
79  GHz.  It  was  shown  that  the  origin  of  the  chirp  is  in 
the  dye-laser  oscillator,  rather  than  in  the  amplifier 
chain. 
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APPENDIX  B 


FOUR  PHOTON  RESONANT  COHERENT  INTERACTION 

We  consider  an  atomic  system  excited  by  a  coherent 
electromagnetic  pulse  which  is  in  condition  of  four  photon 
resonance  between  the  ground  state  |1>  and  a  state  |  2>.  The  6 1 S0 
-  6-*-D2  transition  of  mercury  will  be  used  as  a  numerical  example. 
The  calculations  for  the  6s  -  7s  transition  are  slightly  simpler, 
since  the  energy  of  a  single  photon  is  not  sufficient  to  make 
transitions  from  the  upper  level  (7s)  to  the  continuum. 


Let  us  designate  by  {|^>}  the  set  of  states  that  are  coupled 
to  |l>  or  |  2  >  by  dipole  moments.  It  is  postulated  that  no  other 
transition  than  |1>  -  |2>  is  resonant  with  the  applied  field.  The 
adiabatic  approximation  can  be  used  to  solve  the  equations  of 
motion  for  the  elements  of  the  density  matrix  connecting  a  state 
|^>  to  the  level  |1>  of  \2>  The  case  °f  the  state 
\i>  in  the  continuum  has  to  be  treated  differently,  since  the 
transition  |2>  -  1 1>  can  then  be  considered  to  be  exactly 
resonant.  A  perturbation  treatment  is  used  in  that  case, 
justified  by  the  weakness  of  the  bound-free  matrix  elements.  It 
should  be  noted  however  that  the  latter  approximation  has  never 
been  tested  in  the  femtosecond  range,  and  that  experiments  will 
be  made  to  test  this  particular  aspect  of  the  model. 
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The  linearly  polarized  electromagnetic  pulse  is  described  by 
the  electric  field  vector: 

E=E(t)eitJt+c.c.  (1) 

where  E(t)  is  the  slowly  varying  complex  amplitude  of  the  field. 
The  equations  of  motion  for  the  density  matrix: 

I?  *  hi  <2> 

can  be  solved  by  adiabatic  approximation  for  the  density  matrix 
elements  involving  the  levels  |/>,  essentially  reducing  the 
description  of  the  atomic  system  to  that  of  a  two-level  system: 
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In  Eqs  (3)  -  (5),  j  are  the  slowly  varying  amplitudes  of  the 
density  matrix  elements.  The  transition  matrix  element 
responsible  for  four  photon  absorption  is: 


r,2  =  -  **lf  ^fk  l  _ 

„  <  3w  *  ut\  )  (  2u>  ~  ^kl >  ( w  ~ 

f4 

In  the  case  of  the  6s  -  6d  transition  of  mercury: 
r,,  »  -  2.03  10"162(C-m)4s3 


The  detuning  is  defined  as: 


AtJ  =4 to  -  0>21  -  iw 21 


where  we  include  a  quadratic  Stark  shift  of  the  4-photon 
resonance : 

<5^21  =  )  -  a'(to)]  |e|2  (8) 

In  Eq.  (8),  ce'j_(to)  is  the  real  part  of  the  polarizability  of  the 
state  |i>  at  the  frequency  jj.  The  complex  polarizability  of  |1> 
or  1 2  >  is  given  by: 
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The  proxirriity  of  the  6p  level  resonantly  enhances  the  second  term 
of  the  polarizability  for  level  |  2  >  in  the  case  of  the  6s  -  6d 
transition.  The  values  for  that  transition  are: 

C^U’)  =  4.7  10~40  [  MKS  ]  (1C 

CC2  (to)  =  5.3  10“39  -  i  7.2  lO"40  [MKS]  (11) 

The  imaginary  part  of  the  polarizability  of  state  \ 2>  accounts 
for  the  ionization  rate  from  the  upper  state  (in  Eqs.  (3)  and 
(  4)  )  : 


Y2  =  -2-^2  —  |  E  |  2  (12) 

The  calculation  of  the  time  dependent  macroscopic 
polarizations  at  the  various  frequencies  is  straightforward  but 
tedious,  because  of  the  large  number  of  terms  involved.  For 
instance,  in  the  calculation  of  the  polarization  at  to  ,  one  has 
to  include  intensity  dependent  one-photon  polarizabilities  of  the 
form  )E|2<722E'*  IeI2<7he;  Ie*!3<712'  in  addition  to  the  term  in 
ai  <*iiE  that  appeared  already  in  the  two-photon  calculations  3  • 

In  the  case  of  the  6s  -  6d  transition  of  mercury,  the 
polarization  at  u>  reduces  to: 
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The  polarizations  have^  been  calculated  at  the  various 
frequencies,  and  inserted  in  Maxwell's  equations,  to  compute 
numerically  the  propagation  of  the  various  fields  and  the 
conversion  efficiencies.  The  case  of  the  6s  -  6d  transition  is 
particularly  complex  because  of  the  large  Stark  shift.  In 
addition,  ionization  from  the  upper  level  tends  to  wash  out  all 
coherent  interaction  effects,  except  for  the  shortest  pulses,  for 
which  the  validity  of  our  model  has  to  be  tested. 


The  6s  -  7s  transition  appears  to  be  an  ideal  candidate  for 
efficient  fifth  harmonic  generation.  The  proximity  of  a  9p  level 
near  the  fifth  harmonic  frequency  accounts  for  a  resonant 
enhancement  of  the  fifth  harmonic  generation.  Low  loss 
propagation  of  properly  phased  pairs  of  pulses  is  made  possible 
by  the  absence  of  single  photon  photoionization  from  the  upper 
level  (7s)  and  the  reduced  stark  shift  of  the  resonant 
t  ransit ion . 


APPEND! X  C 


DESIGN  OF  A  50  PSBC  LINEAR  LASER,  AND  ANALYSIS  OP  THE 
COLLIDING  POLSE  EFFECT. 

Jean-Claude  Diels,  Nooshin  Jamasbi 
Center  for  Applied  Quantum  Electronics,  P.O.  Box  5368 
North  Texas  state  University, 

Denton,  Texas  76203.  Ph  (817)  565  3261. 

and 

Laurent  Sarger 

Centre  de  Physique  Moleculaire  Optique  et  Herzienne, 

Universite  de  Bordeaux  I, 

351  cours  de  la  Liberation,  33405  Talence,  France. 

An  antiresonant  (AR)  ring  has  been  proposed  by  siegmari  ^  to  obtain 
standing  wave  saturation  in  a  linear  mode-locked  dye  laser.  This  laser 
configuration  is  equivalent  to  that  of  a  mode-locked  ring  laser  if  the 
saturable  absorber  is  located  at  equal  optical  distance  from  the  beam 
splitter  (Fig.  1).  Two  technical  problems  seem  to  make  this 
configuration  impractical: 

a)  The  requirement  of  tight  focusing  in  the  absorber  2/  3  implies 
a  large  beam  size  through  the  beam  splitter,  hence  a  large 
angle  of  the  AR  triangle,  and  a  large  astigmatism  impeding  the 
tight  focusing. 

b)  it  seems  impractical  to  fulfill  the  "colliding  pulse"  condition 
for  femtosecond  pulses,  since  it  requires  to  maintain  equality 
of  the  two  arms  of  the  AR  ring  with  an  accuracy  of  10  microns. 

We  demonstrate  that  the  astigmatism  of  the  AR  ring  (even  for 
angles  as  large  as  16  °)  can  be  compensated  by  the  geometry  of  the 
cavity  (the  astigmatism  introduced  by  the  other  arm  of  the  cavity  can 
be  used  to  compensate  that  of  the  antiresonant  ring).  Furthermore,  this 
cavity  configuration  enables  us  to  test  for  the  first  time  experi¬ 
mentally  the  "colliding  pulse"  condition,  independently  of  any  other 
parameter.  It  is  shown  that  the  accuracy,  requirement  for  the  equality 
of  the  two  arms  of  the  AR  ring  is  much  less  stringent  than  anticipated. 
The  operation  of  the  laser  is  not  affected  by  a  departure  from  the 
"ideal"  central  position  of  less  than  400  microns! 


The  laser  cavity  is  simple  enough  to  make  a  complete  analytical 
ABCD  matrix  analysis  possible  (using  the  algebraic  formula  manipulation 
computer  language  "MACSYMA").  For  any  value  of  the  angle  0,  we  compute 
the  values  of  mirror  spacing  D  and  2z  that  satisfy  the  stability  crite¬ 
ria  in  the  plane  of  the  cavity  and  in  the  orthogonal  plane.  The  resul¬ 
ting  stability  diagram  shows  "forbidden"  ranges  of  angle  0  for  ’which  no 
stable  cavity  can  be  designed  (for  instance,  13.5°  <  0  <  17°). 
computation  of  the  beam  waists  in  the  absorber  indicate  a  large  spot 
size  (50  microns  and  more)  for  the  first  stability  region.  Astigmatism 
compensation  is  optimized  in  the  small  isolated  stability  range 
corresponding  to  114  <  D  <  117mm,  for  21°  <  9  <  23°,  where  round  spots 
a  few  microns  in  diameter  are  obtained.  Because  of  the  restricted 
dimensions  of  that  stability  region,  the  parameters  of  the  laser  cavity 
are  accurately  defined.  For  a  given  angle  0  (21°  <  0  <  23°),  the  beam 
waist  (in  the  absorber)  is  a  monotonous  function  of  the  mirror  spacing 
D,  making  it  possible  to  study  the  operation  of  the  laser  a£  a  function 
of  absorber  saturation. 

All  the  parameter  (angles,  mirror  distances,  positions  of  the 
images  of  the  fluorescence  spots)  of  the  laser  having  been  predeter¬ 
mined,  the  laser  alignment  follows  a  systematic  and  straightforward 
procedure.  Using  a  0.7  mm  thick  beam  splitter,  a  very  small  amount  of 
positive  dispersion  is  needed  to  compensate  the  small  phase  modulation 
in  the  strongly  saturated  absorber  dye,  resulting  in  a  stable  mode- 
locked  train  of  pulses  as  short  as  50  fsec. 

The  "colliding  pulse"  effect  can  be  investigated  by  changing  the 
relative  length  of  the  two  arms  of  the  antiresonant  ring,  as  sketched 
in  Fig.  2.  As  the  beam  splitter  is  moved  away  from  the  position  corres¬ 
ponding  to  the  symmetric  ring,  the  threshold  pump  power  for  short  pulse 
operation  remains  constant  in  a  800  micron  wide  region.  It  is  only  for 
larger  departures  from  the  "equal  arms"  condition  that  this  threshold 
power  is  seen  to  increase  from  1.7  to  3.5  W.  If  the  pump  power  is 
maintained  at  3.5  W,  none  of  the  following  parameters  of  the  laser  is 
affected  by  departing  from  the  symmetric  ring:  the  pulse  output  power, 
the  intensity  autocorrelation,  and  the  losses  of  the  antiresonant  ring. 
Below  the  "threshold  power  for  short  pulse  operation",  the  laser  output 
consists  in  a  train  of  long  (psec)  bursts  of  noise.  It  can  be  concluded 
that  the  "colliding  pulse  operation"  is  important  in  reducing  the 
threshold  for  stable  mode-locked  operation,  and  in  the  formation  stage 
of  the  steady  state  pulse  from  noise. 
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DESIGN  OF  A  50  PSBC  LINEAR  LASER,  AND  ANALYSIS  OP  THE 
COLLIDING  PULSE  EFFECT. 

SUPPORTING  MATERIAL. 

Analysis  of  existing  femtosecond  ring  lasers  led  to  a  fairly 
complete  understanding  of  the  mechanisms  of  short  pulse  generation  and 
compression  in  these  lasers  1,  2.  Accurate  diagnostic  methods  that  have 
been  developed  confirmed  details  {pulse  shape  and  phase  modulation)  of 
the  theory  3.  it  should  thus  be  possible  to  design  a  laser  based  on 
that  theory.  Our  main  requirements  arc: 

-  a  minimum  number  of  incracavity  elements; 

-  a  cavity  simple  enough  to  make  a  complete  parametric  study 
possible; 

-  tight  focusing  in  the  absorber; 

-  standing  wave  saturation; 

-  an  adjustable  repetition  rate. 

Consistent  with  these  requirements,  we  chose  a  linear  laser  cavity 
terminated  by  an  antiresonant  ring  4  (pig.  i).  imposing  a  spot  size  of 
a  few  microns  in  the  absorber  results  in  a  beam  size  of  4  mm  on  the 
focusing  mirrors,  therefore  a  minimum  angle  of  incidence  on  these 
mirrors,  hence  a  large  astigmatism.  The  latter  effect  appears  to 
conflict  with  the  tight  focusing  requirement.  Analytical  calculations 
of  the  empty  geometrical  cavity  show  that  the  astigmatism  at  the 
absorber  jet  can  be  compensated  by  that  at  the  amplifier  arm.  A  stabi¬ 
lity  diagram  is  plotted  in  Fig.  3,  for  radii  of  curvature  of  7.5  cm  at 
the  amplifier  jet,  3  cm  at  the  absorber  jet.  No  stable  cavity  can  be 
designed  for  13.5°  <  9  <  17  o.  While  the  first  stability  region  is  very 
large,  it  is  only  at  the  edge  of  this  region  that  a  minimum  beam  waist 
(in  the  absorber)  of  50  microns  can  be  calculated.  In  the  very  small 
stability  region  corresponding  to  114  <  D  <  117  mm  and  21  °  <  ©  <  23  °, 
the  beam  waists  are  circular,  with  a  minimum  diameter  of  a  4  microns. 
The  calculation  of  the  cavity  is  complemented  by  calculations  of  the 
size  and  shape  of  the  beam  in  the  amplifying  jet.  Because  of  the  large 
astigmatism  in  that  arm  of  the  laser,  the  pumping  geometry  (focal 
distance  of  the  argon  laser  focusing  mirror  and  angle  of  incidence)  is 
calculated  to  match  the  size  and  shape  of  the  dye  laser  beam  in  the 
amplifier  jet. 


The  (empty  cavity)  laser  has  a  threshold  of  400  mW  (2%  output 


coupling).  With  the  saturable  absorber,  the  laser  exhibits  the  stable 
mode-locking  behavior  typical  of  ring  lasers.  The  pulse  duration  has  a 


sharp  minimum  of  50  fsec  for  an  optimal  amount  of  intracavity  glass. 
However,  iterative  fitting  4  of  the  pulse  spectrum  (fig.  4a),  the 


interferometric  autocorrelation  (fig.  4b)  and  the  intensity  auto¬ 
correlation  (Fig.  4c)  lead  to  an  asymmetric  pulse  shape  with  a  residual 
nonlinear  phase  modulation  (in  contrast  to  the  measurements  with  the 
ring  laser).  It  is  therefore  expected  that  even  shorter  pulses  could  be 
generated  with  this  laser  design,  with  a  better  match  of  the 


intracavity  dispersion  and  phase  modulation. 


Since  we  observe  that  the  two  arms  of  the  AR  ring  can  be  out  of 
balance  by  as  much  as  0.8  mm,  one  may  question  the  need  for  a 
"colliding  pulse"  configuration.  For  the  case  of  astigmatism  compensa- 


tion  at  the  absorber,  the  folded  cavity  obtained  by  removing  the  beam 
splitter  (and  adjusting  z)  is  geometry  identical  to  the  AR  ring  cavity. 
Identical  intensity  autocorrelation  traces  are  obtained  with  that 
cavity,  which  is  essentially  the  same  as  that  used  over  a  decade  ago  5. 
However,  the  stability  is  reduced,  and  the  threshold  a  factor  two 
higher.  r  I  ■  ■  I  ;  1 1-  .  ! 


In  conclusion,  we  have  successfully  designed  a  laser  cavity 
for  which  a  straightforward,  simple  alignment  procedure  replaces 
the  skills  and  artful  expertise  normally  required  for  fsec  lasers. 
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The  use  of  pulse  sequences  offers  a  solution  to  the  conflicting  requirements  of  spectral  selectivity 
and  high  peak  intensities  needed  to  study  multiphoton  processes.  In  addition,  it  is  shown  that 
random — rather  than  regular — pulse  sequences  can  provide  unique  information  on  the  influence  of 
coherence  on  a  resonant  multiphoton  process,  because  the  coherence  can  be  adjusted  independently 
>f  the  resonance  condition. 


I.  INTRODUCTION 

Attention  has  been  drawn  in  recent  years  to  the  possible 
use  of  pulse  sequences  to  enhance  quantum  beat  phenome¬ 
na  tn  molecules,1  or  to  selectively  enhance  weak  spectral 
lines.2  Phase  shifts  in  the  light  field  for  a  single  pair  of 
pulses  can  be  used  to  distinguish  between  single  and  mul¬ 
tiphoton  processes.3  Statistically  distributed  phase 
switching  on  a  cw  field  can  give  it  the  spectral  and  driv¬ 
ing  characteristics  of  a  line-broadened  source.4  Here  we 
discuss  an  application  of  pulse  sequences  in  the  area  of 
laser-driven  reactions.  We  choose  as  an  example  multi¬ 
photon  ionization  of  atoms.  However,  the  approach  used 
here  applies  also  to  multiphoton  dissociation  of 
molecules — which  has  been  the  subject  of  extensive 
theoretical  investigation5 — or  to  any  other  process  depen¬ 
dent  on  the  absorption  of  several  photons. 

Two  types  of  sequential  excitation  will  be  considered: 
deterministic  and  statistical.  The  motivation  for  using  a 
sequence  of  identical  pulses  with  a  well-defined  phase  re¬ 
lationship  is  to  combine  the  spectral  selectivity  offered  by 
a  regular  pulse  train  with  the  high  peak  intensities  of  pi¬ 
cosecond  pulses,  which  are  needed  to  excite  higher-order 
processes. 

A  dimension  of  nonlinear  mechanisms  that  is  subject  to 
experimental  test  is  the  coherence  or  lack  of  coherence  in 
the  process.  But,  as  we  discuss  below,  the  single-pulse  ex¬ 
periments  that  have  been  the  mainstay  of  this  area  provide 
a  fairly  limited  test  of  coherence.  Multimode  versus 
single-mode  operation  lias  been  used"''  to  vary  the  "coher¬ 
ence”  of  the  laser  source.  This  procedure  has  the  disad¬ 
vantage  of  varying  the  bandwidth — hence  the  resonance 
condition — simultaneously  with  the  pulse  coherence.  We 
will  show  how  a  sequence  of  pulses  of  controlled  phase 
can  be  used  to  vary  the  source  coherence  without  affecting 
its  central  frequency  or  its  spectral  width. 

By  choosing  the  delays  and  phase  between  pulses  to  fit 
a  chosen  statistical  distribution  we  can  study  how  a  multi¬ 
photon  process  is  affected  by  those  statistics.  Such  an  ap¬ 
proach,  of  averaging  many  measurements  with  a  predeter- 
•'■l’\s  ''lit  nlviws  Hi'  lit  I  "!>  '.'II  I'hoton 


statistic),  will  be  shown  to  have  several  advantages  over 
other  methods  of  studying  the  influence  of  source  coher¬ 
ence  on  a  multiphoton  resonant  process.  As  an  illustra¬ 
tive  example  for  the  latter  we  have  chosen  multiphoton 
ionization  in  a  metal  vapor  in  the  presence  of  a  two- 
photon  resonance. 

II.  THEORETICAL  MODEL 

We  consider  an  atomic  system  in  conditions  of  two- 
photon  resonance  with  electromagnetic  radiation  at  w, 
which  can  be  ionized  with  a  single-photon  transition  from 
the  upper  state  (Fig.  1).  The  evolution  equations  for  the 
components  of  the  density  matrix  for  an  exciting  electric 
field  of  slowly  varying  amplitude  iff  can  be  written8 


icrl2~  ~ 


17’  i 

iAw+  —  +  ~  / cr , 2  —  (<r:2  —  au)tf * 

Ti  2 


(1) 


dn  —  On=4Re{ff*:(/t7|2)]  —  T’O:;  , 


(2) 


CT j  ,  +(7;;=  —  ,  <3) 

where  is  the  phase  relaxation  time  and  Yi  =  2y6  '~  the 
photoionization  rate  from  the  upper  level  2,  and 
&<o  =  2a>  —  (Ji;  —  is  the  frequency  detuning  from  reso¬ 
nance,  including  the  Stark  shift  bh)\z  —  a6  ~ .  With  the 
substitution  iO’,2~Q<  (o^—  o,|1)/2=  IF,  and  E~  =  2<5;, 
and  neglecting  relaxation  and  ionization,  the  equations 
can  be  rewritten 


Q  =  -iAioQ  -  WE1  .  (4) 

W'-Re<£:*£>>  (5) 

with  the  initial  condition  IFir  =0)=  —  1. 

Equation  (41  can  be  formally  integrated,  yielding 

'■*“"  =  -  ['  WE'e'^'dt'  .  ibi 
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After  a  pulse  has  passed,  {>(/ »r)  =  £»/.(  11),  where 
£:,( lit  is  the  Fourier  transform  of  £:(f).  From  the  con- 
sevation  re’ation  Q  ~  ‘  14'--!, 


2r)S 


J.-C  DIELS  AND  J.  STONE 


31 


Q  -^l-H  :=_2H’0AIK=2AH'  ,  (7) 

where  AH'=  H'(f  =  x  j-~  IK(r=0).  Combining  Eqs.  <6' 
and  1 7)  gives  the  energy  absorbed  by  the  two-!evel  system 

2  AH’  =  —  £_V(A«>;2.  'SI 

The  number  of  ions  produced  is  simply  the  number  of 
"lost  atoms"  1  —  =  -r. )  —  a{^t  =  x )  in  the  system  of 

Eqs.  (1) — 3).  The  number  of  ions  can  also  be  estimated  by 
integrating  the  ionization  rate 

-V.,,,,s  =  f  y  E  zcr:idt  .  (9) 

For  small  population  transfers,  <7::  —  AH’  and 

•V.un,=  /r'£i:A»*(rW  r.  ilO> 

Replacing  AIK  by  its  average  value  W{  x  )/2, 

Mns  —  7  Y  i  E->py  At >)  ■  ~9  ,  (11) 


<•  *  , 

where  9  =  J  \E  \  ~dt  is  the  "pulse  area"  or  tipping  an¬ 
gle  of  the  pseudopolarization  vector  in  the  vector  model 
of  the  two-photon  resonant  interaction.’10 

The  expressions  above  hold  in  the  approximation  of 
weak  excitation  ( 6  small),  no  Stark  shift  (or  any  Stark 
shift  but  square  pulse  excitation).  We  see  that  both  the 
absorption  and  the  photoionization  reflect  the  excitation 
spectrum.  The  photoionization  spectrum  is  simply  a  frac¬ 
tion  (proportional  to  the  pulse  area)  of  the  Fourier  spec¬ 
trum  of  the  Rabi  frequency  (which  is,  in  this  case,  propor¬ 
tional  to  the  square  of  the  exciting  electric  field). 

III.  COHERENCE  AND  SPECTRAL 
CHARACTERISTICS 

We  will  consider  now  more  specifically  the  spectral 
characteristics  (i.e.,  number  of  ions  versus  frequency)  of 
an  excitation  mode  of  a  sequence  of  pulses.  In  general, 
even  if  each  individual  pulse  is  pure  in  frequency  (real 
amplitude  function)  the  sequence  of  various  pulses  will  be 
represented  by  a  complex  function  E  (because  of  the  rela¬ 
tive  phase  between  pulses).  A  general  property  of  the  sys¬ 
tem  of  Eqs.  (4)  and  (5)  is  that  H’ — which  determines  in 
first  approximation  the  photoionization — remains  un¬ 
changed  under  a  transformation  £—»£*,  Aoi— ►  —  Aa» 
(indeed,  under  this  transformation  Eq.  (4)  is  transformed 
in  its  complex  conjugate  and  Eq.  (5)  remains  unchanged]. 
This  symmetry  implies  that  the  photoionization  spectrum 
of  a  single  unchirped  pulse  (£  real)  is  symmetric.  The 
symmetry  is  broken,  however,  in  the  case  of  a  sequence  of 
pulses  with  a  phase  difference  A <b  different  from  0  to  rr, 


O 
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FIG.  2.  Plot  of  (he  spectral  intensity  for  a  sequence  of  ten  identical  square  pulses,  as  defined  in  Eq.  05). 
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that  of  the  excitation  pulses.  In  fact,  the  photoionization 
curves  show  all  the  characteristic  features  of  a  Fourier 
spectrum:  the  overall  width  is  related  to  1/r,  an  increase 
of  the  phase  difference  between  pulses  shifts  correspond¬ 
ingly  the  photoionization  spectrum,  and  the  width  of  the 
individual  peaks  in  the  photoionization  spectrum  de¬ 
creases  with  incre..'ing  delay  between  pulses. 

The  Stark  shift  introduces  an  additional  complexity  by 
essentially  creating  a  different  resonance  condition  during 
and  between  pulses.  This  effect  is  demonstrated  by  the 
simple  example  of  a  sequence  of  two  square  pulses  shown 
in  Fig.  4.  Increments  of  delay  between  the  pulses  by  mul¬ 
tiples  of  2T/iStark  shift)  produce  the  same  number  of  ions 
at  resonance  and  at  a  detuning  from  resonance  equal  to 
the  Stark  shift.  Departure  from  t he  simple  approxima¬ 
tions  [Eqs.  (6i— ilOi]  will  be  particularly  significant  for 
large-area  pulses.  "Photoionization  spectroscopy”  with  a 
sequence  of  large-area  pulses  will  enable  a  direct  measure¬ 
ment  of  the  photoionization  rate  from  the  upper  state 
through  the  width  of  ihe  individual  peaks.  Indeed,  as 
shi-  vn  in  F:g.  5,  ihe  width  of  each  peak  in  die  photoioni- 
zation  spectrum  is  seer  to  increase  with  the  transition  r  ite 
to  the  continuum.  It  is  interesting  to  note  Fig.  5i  that  the 
peak  ionization  is  swjtier  for  the  higher  ionization  rate. 
The  reduction  in  peak  ionization  rate*  is  due  to  the  in¬ 
creased  dephasing  associated  with  the  transition  rate  to 
the  continuum,  resulting  in  a  decrease  in  coherent  excita- 
'ion  of  the  upper  level. 


PULSE  "AREA";3 


FIG.  ?.  Single  peak  of  the  ionization  spe*ctrum  G-  ions  vs  an¬ 
gular  frequency  in  nsec")  for  two  different  ionization  rates. 
The  exciting  field  consists  of  a  sequence  of  ten  Gaussian  pulses, 
of  0.01  nsec  duration  TWUM*.  and  area  O-  j  E:dt  -  3.  The 
successive  pubes  are  0  4  oul  of  phase.  The  medium  parameters 
are  the  same  as  in  the  preceding  figures,  except  for  the  ioniza¬ 
tion  rates  yE:  being  considered  here.  The  scale  to  the  right  cor¬ 
responds  to  the  ionization  rate  of  0.07 /•.'•  i  -  -  -  I,  while  ihe 
scale  to  the  left  pertains  to  the  higher  ionization  rate  of  0  44£; 


V.  IONIZATION  BY  A  RANDOM  PULSE  SEQUENCE 

An  incoherent  field  can  be  described  by  a  density  ma¬ 
trix  corresponding  to  an  incoherent  superposition  of 
coherent  states."  We  choose  to  use  for  this  superposition 
a  series  of  identical  pulses  with  a  random  distribution  of 
delays  and  phases. 

Under  certain  conditions  a  large  number  of  experiments 
performed  with  a  few  pulses,  for  which  the  relative  delays 
are  distributed  according,  for  instance,  to  a  Poisson  statis¬ 
tics,  will  have  the  same  average  result  as  the  same  experi¬ 
ment  perform- J  once  with  a  very  large  number  of  pul  .es, 
distributed  according  to  the  same  statistics.  The  average 
time  between  pulses  has  to  be  shorter  than  the  coherence 
time  of  the  particular  resonance  under  investigation.  We 
demonstrate  next  that,  if  the  pulses  are  randomly  distri¬ 
buted,  the  time  evolution  of  the  system  driven  by  a  se¬ 
quence  of  identical  pulses  is  proportional  to  the  number  of 
pulses  and  to  the  total  fluence  delivered,  as  expected  for 
incoherent  excitation.  If,  instead,  the  pulses  are  applied  in 
phase  and.  for  instance,  in  the  condition  of  two-photon 
resonance,  the  absorption  by  the  system  is  proportional  to 
the  square  of  the  number  of  pulses.  For  demonstration 
purposes  let  us  consider  near  two-photon  resonance 
(£:»3o)  and  weak  pulses  such  that  the  approximation 
1F~IF)  can  be  made  in  Eq.  <4).  Each  of  the  pulses  of 
“area”  30  =  J  '  E  dt  produces  an  increase  in  polarization 

O  =  IFoA0exp(i(M  .  i  16) 

If  the  pulses  are  all  in  phase  pulses  equally  spaced  and  in 
phase),  30  is  a  constant  angle,  and  after  .V  pulses, 
Q  =  .Y30  and  the  energy  absorbed  IV  is  proportional  to 
the  square  of  the  number  of  pulses: 


If  the  distribution  of  pulse  delays  is  random  1  Surj  «  I ), 
then  the  phase  in  Eq.  ( 1)  is  a  random  variable  and  the  vec¬ 
tor  Q  makes  a  random  walk  of  .V  steps  30  in  the  complex 
plane,  starting  from  the  origin.  The  resultant  polarization 
is  Q=\ 30.  and  the  energy  absorbed  is  proportional  to 
the  number  of  pulses: 


It  should  be  noted  that  by  using  a  sequence  of  pulses  we 
can  compare  coherent  [as  in  Eq.  17>)  versus  incoherent 
[as  in  Eq.  (IS)]  excitation,  under  the  same  conditions  of  to¬ 
tal  fluence  and  peak  intensity.  As  expected,  it  is  only  in 
the  case  of  incoherent  excitation  that  the  evolution  of  the 
system  can  be  simply  proportional  to  the  total  fluence 
delivered.  The  intensity  dependence  of  the  multiphoton 
process  is  not  invoked  in  comparing  the  incoherent  versus 
coherent  ease.  The  above  considerations  still  hold  in  the 
more  general  case  of  a  n-photon  transition.  Ihe  only 
difference  is  a  generalized  Rabi  frequency  in  £”.  Equa- 

f 

E"dt 

In  the  approximation  En  »3ro,  and  no  intermediate  level 
resonance,  coherent  excitation  is  characterized  by  an  ab¬ 
sorbed  energy  proportional  to  the  square  of  the  number  of 
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pulses,  while  the  absorption  is  directly  proportional  to  the 
number  of  pulses  in  the  incoherent  case. 

In  order  to  study  the  influence  of  source  coherence  on 
the  photoionization  spectrum  of  lithium  we  calculated  the 
number  of  ions  as  a  function  of  pulse  wavelength  for  a 
coherent  and  an  incoherent  sequence  of  ten  pulses  t  Fig.  6). 
The  "incoherent"  curse  (solid  line.  Fig.  6)  corresponds  to 
an  average  of  500  pulse  sequences,  each  with  a  different 
set  of  relative  delays  taken  from  a  Poisson  distribution. 
The  relative  phases  are  uniform  random  numbers  between 
0  and  2~.  The  total  area  in  each  sequence  of  pulses  is  1. 
As  can  be  seen  from  Fig.  6.  the  photoionization  from  a 
coherent  sequence  i ton  pulses  equally  spaced,  relative 
phase  ~  2,  dashed  line  in  Fig.  tv  can  be  much  larger  than 
that  of  the  incoherent  excitation.  This  result  may  appear 
to  be  in  contradiction  with  experimental  results  of 
Lecompte  cl  a!.."  on  the  influence  of  coherence  on  multi¬ 
photon  excitation.  In  the  case  of  nonresor.ant  .V  photon 
ionization  of  atoms  reported  in  Ref.  b,  the  ionization 
probability  is  proportional  to  the  .Vth  power  of  the  laser 
intensity  and  therefore  scry  sensitive  to  the  transient  in¬ 
tensity  peaks  that  occur  in  a  multimode  laser  pulse.  The 
dependence  of  multiphoton  ionization  on  laser  coherence 
was  studied  by  varying  the  number  of  oscillating  laser 
cavity  modes  between  me  and  one  hundred.  It  has  been 
•est'.mated"  that,  tn  that  case,  a  20-mode  laser  is  a  good 
approximation  of  an  incoherent  source.  The  situation  is 
more  complex  in  the  case  of  resonant  multiphoton  ioniza¬ 
tion  as  we  are  considering  here.  Indeed,  the  number  of  os- 
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FIG.  6.  Ionization  ■. s  angular  frequency  msec  1  for  a  se¬ 
quence  of  ten  pulses.  Each  pulse  is  Gaussian  shaped  of  area  0. 1 . 
The  solid  line  corresponds  to  a  random  pulse  sequence,  the 
dashed  line  to  a  sequence  of  equally  spaced  '0.09  nseei  pulses, 
each  90'  out  of  phase  with  the  former.  F.acli  pulse  of  the  se¬ 
quence  is  Gaussian,  with  a  duration  of  0.01  nsec  FWIIM1.  The 

statistical  result  1 - i  is  obtained  by  taking  the  average  of  500 

"excitations."  each  with  a  sequence  of  pulses  distributed  'in  rela¬ 
tive  delay1  according  lo  a  Poisson  distribution.  The  relative 
pluses  between  pulses  are  taken  lo  be  random  numbers  between 
0  and  2-  The  medium  parameters  are  the  same  as  for  Fig.  2. 
The  average  delay  between  pulses  is  0.09  nsec. 


cillating  modes  is  varied  by  changing  the  bandwidth  of 
the  cavity.  Therefore,  the  parameter  coherence  cannot  be 
controlled  independently  of  the  resonance  condition,  since 
the  spectral  width  of  the  laser  source  increases  with  de¬ 
creasing  coherence.  In  the  approach  used  here,  however, 
the  laser  coherence  is  modified  at  constant  bandwidth. 
since  the  spectrum  of  the  pulse  sequence  is  contained 
within  the  spectrum  of  a  single  pulse. 

It  is  interesting  to  note  that  in  a  recent  experimental 
study  of  four-photon  three-photon  resonant)  ionization  of 
caesium  atoms,  Lompre  et  alP  found  also  an  increase  in 
the  width  of  the  spectral  dependence  of  the  ionization  for 
incoherent  versus  coherent  excitation.  The  hyperfine 
structure  of  the  resonance  is  still  observed  in  coherent  ex¬ 
citation,  while  it  is  not  in  the  incoherent  case  (Ref.  7,  Fig. 
8).  However,  in  that  particular  experiment7  the  coherence 
was  varied  by  changing  the  numbe'r  of  modes  of  the  laser 
and  comparing  data  for  the  same  average  intensity. 
Therefore,  temporal  fluctuations  make  the  peak  intensity 
highest  in  the  incoherent  case,  resulting  in  a  higher  ioniza¬ 
tion  in  the  incoherent  case  than  in  the  coherent  case.  In 
an  experiment  in  which  the  coherence  is  varied  by  varying 
the  temporal  distribution  of  identical  pulses,  as  suggested 
here,  the  peak  intensity  as  well  as  the  spectral  bandwidth 
remain  constant,  while  the  coherence  is  modified  within 
the  pulse  spectrum. 

Photon  statistics.  Let  us  consider  a  limiting  case  where 
each  pulse  of  the  sequence  contains  only  one  photon.  The 
"pulse  sequence  statistics’*  is  then  nothing  more  than  the 
photon  statistics  of  the  radiation  field.  The  number  of 
pulses  V  is  now  the  number  of  photons  in  the  excitation, 
each  of  these  photons  fiu  driving  coherently  an  n-photon 
resonant  system  i.e.,  a  two-level  system  of  energy  n  fuj). 
For  instance,  in  the  simple  argument  at  the  beginning  of 
this  section,  the  elementary  increment  of  the  n-photon 
off-diagonal  clement  becomes 

A0  =  ['/  Endt=T]if,(o)nr~  ,  (19) 

7/  being  a  proportionality  coefficient. 

It  should  be  noted  that  the  photon  energy  does  not  ap¬ 
pear  explicitly  in  any  calculations  carried  out  in  the  slow¬ 
ly  varying  envelope  approximation  (SVEA).  The  statisti¬ 
cal  sequence  of  identical  macroscopic  coherent  pulses 
(each  of  them  being  a  coherent  superposition  of  Af  pho¬ 
tons)  interacting  with  the  n-photon  resonant  atoms  is 
identical  (within  the  SVEA)  to  a  "photon  statistics’*  in  a 
fictitious  world  where  Planck’s  constant  is  Af ft.  The  tran¬ 
sition  frequency  at  is  in  the  condition  of  «- 

photon  resonance  with  each  pulse  ( Mfi)e>n~n  (Mfuo). 
Since  Planck’s  constant  appears  only  implicitly  (in  the 
scaling  of  the  generalized  Rabi  frequency  E)  in  the 
theoretical  model,  all  the  equations  of  the  theoretical 
model  remain  unchanged  by  the  "transformation."  The 
elementary  increment  of  the  off-diagonal  element  for  the 
standard  pulse  of  the  train  can  also  be  expressed  in  terms 
of  the  "pseudophoton”  energy  (Affl ii».  in  a  similar  way  as 
in  Eq.  (19): 

A 0=  f  Endt  =)/[(.' .  120) 

This  similitude  offers  an  intuitive  justification  for  the  use 
of  a  pulse  sequence  to  study  the  statistical  response  of  a 
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multiphoton  resonant  system.  We  are  essentially  creating 
a  "pseudophoton  statistics."  This  concept  is  legitimized 
because  Planck's  constant  does  not  appear  explicitly  in 
this  formalism  ti  e.,  the  photon  energy  is  still  small  com¬ 
pared  with  the  total  excitation  energy).  Since  the  model 
of  the  pseudostatistics  concept  is  equivalent  to  the  theoret¬ 
ical  model  used  in  Sec.  III.  it  does  not  bring  any  new  light 
to  the  problem.  However,  it  offers  a  convenient  interpre¬ 
tation  to  the  approach  that  is  being  taken  here,  to  study 
coherence  effects  through  the  use  of  pulse  sequences. 

VI.  experimental  considerations 

In  this  section  we  briefly  discuss  the  experimental  im¬ 
plementation  of  the  excitation  scheme  discussed  above. 
The  use  of  a  Mach  Zehnder  interferometer  to  generate 
two-pulse  sequences  has  been  demonstrated  in  a  study  ot 
two-photon  resonant  coherent  propagation  in  lithium  va¬ 
por.1314  Optical  delays  adjusted  with  interferometric  ac¬ 
curacy  were  used  to  control  the  relative  phase  and  delay 
between  pulses.  A  He-\'e  laser  was  sent  through  the  same 
Mach  Zehnder  interferometer  to  monitor  the  optical  de¬ 
lays  (hence  the  phases)  with  subwavelength  accuracy.  A 
series  of  optical  delay  lines  can  generate  multiple-pulse  se¬ 
quences.  Instead  of  Mach  Zehnder  interferometers  a  cas¬ 
cade  of  .V  optical  delay  lines  as  sketched  in  Fig.  7  can  be 
used  to  generate  a  3.V-pulse  sequence.  For  a  beam-splitter 
reflectivity  of  R  =  68rr,  each  "module"  of  Fig.  7  produces 
a  series  of  three  identical  pulses,  each  carrying  22 %  of  the 
input  pulse  energy. 


FIG.  7.  Sketch  of  an  interferometric  delay  line  to  split  an  in¬ 
put  pulse  into  a  sequence  of  three  identical  pulses.  If  the  beam 
splitters  A  and  B  have,  respectively,  a  reflection  coefficient  of 
0.68  and  0.32,  the  successive  pulses  follow  the  paths  AB,  AL  B. 
and  ABLALB  and  each  carries  22 cc  of  the  original  beam. 


The  shorter  the  wavelength,  the  more  difficult  it  is  to 
control  the  relative  phase  of  pulses  with  optical  delay 
lines.  As  noted  earlier,  the  approach  proposed  here  ap¬ 
plies  also  to  multiphoton  dissociation  of  molecules  involv¬ 
ing,  for  instance,  the  use  of  picosecond  COs  lasers,15  for 
which  coherent  pulse  sequencing  li.e.,  pulse  sequencing 
while  maintaining  the  relative  phases)  would  be  particu¬ 
larly  easy. 
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